NTER MITWIRKUNG VON L. COLLATZ : H.GÖRTLER - J. HEINHOLD . K. 'KLOTTER 
-MARGUERRE . H. NEUBER » L. SCHMETTERER . K. SCHRÖDER - H. SCHUBERT 
‚. TOLLMIEN » H. UNGER UND C. WEBER : HERAUSGEGEBEN VON H. HEINRICH, DRESDEN, 


UND 41 “= SEITE 185—228 HEFT S MAI 1961 


AUS DEM INHALT: 


HUA: UP SAU RS AT Zee 


B. Tanimoto: The Reciprocal Theorem Method 

for Applied Mechanics / D. Rüdiger: Eine geo- 

metrisch nichtlineare Schalentheorie / H. Leip- 
holz: Eine Begründung für die Anwendung der 

Laplace-Transformation auf nichtlineare Differen- 

tialgleichungen 


KELSEIPN EM. FT ST E.EESUENE GEBEN 
B-UC -A..B.:E%S PARSE EHSUUNDGSEES 
EINGEGANGENE BÜCHER 
NACHHER 1 TC EHE 


BEN ERLAG GMBH: - BERLIN 


MM Bd. 41 Nr, 5 S. 185—228 Berlin, Mai 1961 


EEE SEITE IT Eu eur 


INHALT: 
Hauptaufsätze Seite 
B. Tanimoto: The Reciprocal Theorem Method for e 
Applied Mechanics Be ee N; 1 
D. Rüdiger: Eine geometrisch nichtlineare Schalen- % 
HneOmie den a ne east ee Ferien 1 
H. Leipholz: Eine Begründung für die Anwendung der 
BT STransfaßmation auf nichtlineare Differential- 
gleichungen . v .u en enn NE 
Kleine Mitteilungen 
E. Becker: Die Riffelung eines von Wasserwellen über- 
spülten Sandbodens . . »... rer rn. 215 
S. B. Dutt: Stresses in a Twisted Cylindrical Plug of Iso- 
- tropic Material Inserted into a Rigid Medium. . . » » 216 
S. Basuli: Note on the Bending of a Thin Uniformly 
Compressed Circular Plate under an Eccentric Concen- 
redal.oud Ber a re ee ee inte = 217 
S. K. Datta: Flow of Visco-Elastic Maxwell Fluid through 
Tubes under Exponential Pressure Gradient ... . . . 219 
W. Dück: Zur Abschätzung der Fortpflanzung der Daten- E 
fehler bei linearen Gleichungssystemen nach der Formel 
von Wittmeyer-Colladz. .. 2.2 2. nen. 220 
Buchbesprechungen ...... BT 223 
Eingegangene Bücher... ..... 2.2... 226 
Nachrichten „Un. Se Bra 227 


Wir bitten, alleManuskriptsendungen direkt an den Herausgeber, 
Prof. Dr.-Ing. H. Heinrich, Dresden A 27, Friedrich-Hegel-Str. 31, 
zu richten. Zu den Arbeiten, die als Hauptaufsätze bestimmt 
sind, erbitten wir auf gesondertem Blatt eine kurze Zusammen- 
fassung des Inhalts, nach Möglichkeit in deutscher, englischer 
und russischer Sprache; falls die Übersetzungen nicht geliefert 
werden können, ist wenigstens die Angabe spezieller Fachaus- 
drücke in den verschiedenen Sprachen erwünscht. Die Arbeiten 
sollen in klarer Schrift, möglichst mit Schreibmaschine, weit- 
zeilig und einseitig geschrieben sein und die nötigen Hinweise 
für den Setzer betreffend Schrifttypen (z. B. griechisch, Fraktur), 
Sperrungen o. a. enthalten. Zur Beschleunigung des Drucks 
und zur Vermeidung von Satzfehlern empfiehlt es sich, umfang- 
reiche und unübersichtliche Formelausdrücke durch Einführung 
von Abkürzungen zu vermeiden. Bilder sollen als Tusch- 
zeichnungen auf Transparentpapier (Beschriftung mit Bleistift) 
oder als saubere Bleistiftskizzen ausgeführt sein. Bildunter- 
schriften sind am Schluß des Textes anzufügen. Sollte die Arbeit 
bereits an anderer Stelle verbreitet sein (als Dissertation, For- 
schungsbericht, Manuskriptdruck o. ä.), so ist dies auf der ersten 
Textseite in einer Fußnote anzugeben, Für die Zusammenstellung 
der zitiertenLiteratur, amSchluß der Arbeit, bitten wir, sich einer 
einheitlichen Notation entsprechend folgendem Muster zu bedie- 
nen, z.B. W.Schmeidler, Über die Wärmespannungen in ei- 
nem Körper, ZAMM 28 (1948), S. 54-59 oder G. Bürger- 
meister und H.Steup, Stabilitätstheorie I, 1. Aufl., Berlin 1957, 
Akademie-Verlag, S. 142—147. Die Autoren erhalten von den 
Hauptaufsätzen 75, von den Kleinen Mitteilungen 25 Sonder- 
drucke ohne Berechnung, darüber hinaus bis zu 250 Sonder- 
drucke gegen Berechnung. 


Der Verlag behält sich für alle Beiträge das Recht der Verviel- 
fältigung, Verbreitung und Übersetzung vor, 


EBERLE RESTE ET RETTEN ERTTEEN 


SCHRIFTENREIHE 
DER INSTITUTE 
FÜR MATHEMATIK 


BEI DER DEUTSCHEN AKADEMIE 
DER WISSENSCHAFTEN ZU BERLIN 


7 
| 
; 
BE 


von 28 Autoren, 
Prof. Dr. J. NAAS und . Dr. KURT 


1957. VIIL, 317 Seiten — 22 Abbildungen — 9 Tafeln — gr. 8° — 
DM 38,— (vergr 
Lösung des für eindimen- 


Von Dr. KARL BORKMANN und SIEGFRIED OBERLÄN- 
DER 
1955. 99 Seiten — gr. 8° — DM 12,— (vergriffen) 


Ki A.W. EEROERBEN 

bersetzung Russischen 

(Tisscnschaftiiche Redaktion: Preh Dr. J. NAAS 
1956. 79 Seiten — gr. 8° — DM 5,50 


Wirbelstraßen k 


Von Prof. Dr. Bl. DOLAPTSCHIEW 
1957. 28 Seiten — gr. 8° — DM 3,80 


Die Verbiegung konvexer Flächen 
Von A. W. POGORELOW 
(Übersetzung aus dem Russischen) 
Wissenschaftliche Redaktion: Prof. Dr. E. REMBS 
1957. 135 Seiten — 26 Abbildungen — gr. 8° — DM 18,50 


Einfül in die Theori ionärer Zufallsfunkti 
Von A.M. JAGLOM 
rsetzung aus dem Russischen) 
Deutsche Übersetzung unter wissenschaftlicher Redaktion von 
Dr. HERBERT GOERING 2 
1959. VIII, 177 Seiten — 8 Abbildungen — gr. 8° — DM 24,— 


Tabellen von Exponentialfunktionen und -integralen zur 
Anwendung auf Gebieten der Thermodynamik, Halbleiter- 
theorie und ineti 

Von SIEGFRIED OBERLÄNDER 


1959. VIII, 8 Seiten Text — 142 Tabellenseiten — 1 Abbildung — 
gr. 8° — Ganzleinen DM 35,— (vergriffen) 


Konstruktion ganzer, rationaler und reeller Ordinal- 
zahlen und die diskontinuierliche Struktur der trans- 
finiten reellen Zahlenräume 


Von DIETER KLAUA 
1961. 141 Seiten — gr. 8° — DM 28,— 


Qualitative Methoden beim n-Körperproblem der Him- 
melsmechanik 
Von F. CHILMI 
Übersetzung aus dem Russischen; vom Autor durchgesehene und 
ergänzte Übersetzung. In deutscher Sprache herausgegeben von 
J. O0, FLECKENSTEIN 


1961. 116 Seiten — 4 Abbildungen — 1 Tabelle — gr.8° — 
DM 23,— 


Über eine Bewertungstheorie der Algebren und ihre Be- 
deutung für die Arithmetik 

Von HERBERT BENZ 

1961, 149 Seiten — 5 Abbildungen — gr, 8° — DM 29,50 


Fordern Sie bitte unseren Fachkatalog MATHEMATIK an, der 
Ihnen kostenlos geliefert wird. 


Bestellungen durch eine Buchhandlung erbeten 


AKADEMIE-VERLAG . BERLIN 


ERBETEN 


ZEITSCHRIFT FÜR ANGEWANDTE 
MATHEMATIK UND MECHANIK 


INGENIEURWISSENSCHAFTLICHE FORSCHUNGSARBEITEN 
Band 41 Mai Heft 5 


E ; q | ZAMM 41 (1961) Heft 5, Seite 185—197 
The Reciprocal Theorem Method for Applied Mechanics 


By B. TAanımoro 


Das Problem der Durchbiegung von Balken und Platten wird mittels einer Erweiterung des Rezi- 
prozitätssatzes von Maxwell und Betti behandelt. Bei dieser neuen Methode werden simultane 
Gleichungen, wie sie bei der üblichen Lösungsmethode auftreten, vermieden. Im Falle der Platten ist 
die Benutzung der transponierten Matrix besonders wirkungsvoll. 


The bending of beams and plates is treated by the extended use of the reciprocal theorem of Maxwell and 
Betti. This new method expels simultaneous equations to be encountered in the ordinary method of 
analysis. In the case of plates, the use of transposed matrix is especially efficient. 


PaccMmaTpuBaecTA uU3TuUÖaHHe 6AJIOK H ILJIACTUH IPH IOMOIIM IIHPOKOTO UPHMEeHEHUA 00PAaTHoK 
TeopeMbI MakcBeuA u BeTTu. ITOT HOBblÜ MeTON HUCKJIOYAeT ONHOBPEeMEHHbIe YypaBHeHnAn, 
BCTpeyaeMble IIPUH IPMHMeHeHUN OÖBIYHOTO MeTona aHasımaa. B cAIy4yae NANACTUH IPHMeHeHHe 
TPaHCHOHHPOBAHHOÜ MATPHlBI OKA3bIBAeTCH OCOÖEHHO IAPPERTHUBHEIM. 


1. Introduetory Remarks 


The reeiprocal theorem of MAxwELL and BET1T is well known to be a powerful weapon in 
many fields of mechanies. When, for example, applied to the bending of beams, this theorem 
states 


mPıt+m P+-+mPr=mPıtmP2t+t'+mPa for Fig.1i1...(W), 
and accordingly 
means iorBie:27 Mu. a, 2 Rene 
It follows from this that, 
EEE De CHEN FF EN, R) EN ERDE LIEFELIEH 


provided that x and o are respectively the same values. 
Pr % D Peg a 
p 2 © 


ER 7, 7 7=Ffx,9) 


A 


nV 


N2 72:0, ) 
7: 


Fig.1 Fig. 2 Fig. 3 


The present use of the theorem is the extension of equation (3) to the case of Fig. 3, in 
which P denotes a concentrated load applied at point o = x, n, the deflection for the domain 
0<o<x, and n, the deflection for the domain «<e<1. If, be denoted by /(z, o), then n, 
is given by interchanging x and o, that is by f(e, x). 

In virtue of this, simultaneous equations to be encountered in the ordinary method of 
analysis are eliminated, and hence time and labour is much reduced. In the case of plates, the 


adoption of the matrix form in the course of calculation is especially suited. 
13 


Bann Zu no. 
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2. Bending of Beams 


A Beh cases of the bending of beams will be given. 
The first case is the cantilever as illustrated in Fig. 4. Introducing dimensionless quantities 
0 — x/l and x = £/l, the deflection equation for the domain 0 <eo<zis 
d’n, M»+-- Bil 
da. FÜ 23 aa 


EI being flexural rigidity of the beam. Integrating this equation twice with ae to x, and 


00 


taking (n)e-0o = 0 and (2) Bin 0 into account, we at once obtain 


PB 
HT ET 


Then the deflection n, for the alternative domain «<< o<lis given by interchanging the letters 
» and o; that is 


(3% 0? — 0?) (<Se<&M- : 0000 


PB 
2 SE 


Equations (4) and (5) are in accordance with those given in ordinary textbooks [1]. It is to be 
noticed that in this case no continuity conditions at og = x are necessary. 


Be 


BR eo <Dee .. 6). 


Fig. 4 Fig. 5 


The second example is the simple beam (Fig. 5). Defleetion equation for 7, (O<o<x) 
becomes 


den, M Rt 
== = 0 . 
ine 51 gılt=We (<eo<r 
Integrating this equation, with the condition (7) =o = O at the left end A ofthe beam, 
PP 2 
> 6 Er cr).ah-- nr Se er 


where c(x) is a constant depending on the parameter x. Then by the reciprocal theorem, the 
deflection 7, (« <o< 1) is written 


> 783 


Pl 
n=—en; Ü-DEHEOR] 0 Ser 


The continuity conditions for equations (6) and (7) are 


dn (ding 
TERN . (e., 2 
The first of these equations is identically satisfied by equations (6) and (7), while the second gives 


de(x) c(x) “ 


3x7—2n, 
dx x 


Ex) 
do 

the differential equation above can easily be found to be 

cr) = Axr— 2 — a) 1—x)%, 


A being a mere constant. Substituting this solution into equation (7), and considering the condi- 
tion (mn -ı =Oattheright end B, we obtain A =. 


") de(x) 


where TER since c(x) in equation (6) is independent of o. The general solution of 
x 


[0 7 


3 


I En ZZ 
”s 


“ 
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Hence the wanted deflections 7, and n, become respectively 
ur: 
PR (Big.5) a 2a) 
OP 


This result is in accordance with that given in ordinary textbooks [2]. 


The third example is the beam indicated in Fig. 6. This beam is known to be a eätteallg 
indeterminate beam of one degree. The reaction R, at the left end A may easily be found by 
using the Sun ep reciprocal theorem, HR (1), from which we have 


= SU-mM@+n. 


Then the deflection equation becomes 


dene ee Pl Age 
Fre ET ap, x” (2 +) (Sr: 


Integrating this, with the condition (n,). -o = 0, we have 
2,8 


4 2 3 
c(x) being constant. Then by the reciprocal theorem 
PP 
= 2 
Na a5 MORD Hr co)el R<o<). 


The continuity condition at point o = x, equations (8), affords 
dc(x)  c(x) 
ee ee 14 
dx x su 2% 


the general solution of which is 
cd) = Ar +32 —r). 


dr, 
The conditions at the right end B are (n5). -ı = Q and (m —= 0, which afford A = — 3. 
e=1 
Hence the wanted deflections n, and n, become 
PR 
ee ER 2,9) 0 
“pP (ISIN Te 10): 
wi 2 x —(24 & t g 


This result is in accordance with that given in ordinary textbooks [3]. 
The fourth example is the beam ae in Fig. 7. By the ordinary reciprocal theorem or 
otherwise, we have 
Bet -(t+2), M=--Plal—»)%. 
The deflection equation then becomes 


den, MesPl 


DerannrerN ie —- Al+2W)ol P<eo<. 


13* 


+ 
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Integrating this, with the condition that the left end A is clamped, we simply have 
PR 


m BT (0 < e<x, Fig7) ur re 
so that the defleetion of the alternative domain is 
3 ” 
M= 1- Ber —(i+20)#] «<e<1,FigT) .... (2). 


The result obtained is in accordance with that given in ordinary textbooks [4]. 


The last example is the two-spanned continuous beam, in which for convenience equal span 
is assumed, as given in Fig. 8. By the ordinary reciprocal theorem, we first have 


P 
R,= 74-5249), R= 520-9), Re=— rim. 


- 


The deflection equation for 7, becomes 


den, M Pl 
RT = 4gp@r I Te (<e<x). 
Integrating this, with the condition (n,). =o = 0, we have 
BR 


m=—ap MS) + ec) el (<o<»r); 


and accordingly 


PR 
n=— 5 Mer HAr+edr KF<o<)). 
The continuity condition at 0 = x affords 
dc(x) c(x) 
N RR EI EE 
er r x(6—5%), 


the general solution of which is 
ca) =Axr + (12 —5x)x:. 


The condition (ng). - ı = 0 at the intermediate support B gives A = — 7. 
Hence we obtain 
' PR e 
h=aEr 97-5980 Ad) U<o<H, 
PB (Filed) © 2. -for 


n= 457-7 5dor ker] (K<eo<) 


Our last step is to find n,. To do this we first have (Fig. 9) 


tn Pl 
M;=-—E a) True 
and 
d’n P 
R | 12 = — a 
B mal 1 x(5 — 22), 


E 
> 
r 
f 
» 
ur 
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Rz being a part of R,, which plays as the reaction as if AB were a single-spanned beam. Then 
the remaining part of R,, say R%, which is R, — R, — Ri}, is given by 


ker P P 
BR eRr-nh,- er) rl) z nal). 
Hence the deflection equation for 73 becomes 


dp _ . M Pl 


1 " 
dx? ni Drau Ne ee lo 


or, on integration with the conditions Mel ei, PM; 
and (nz). = ı = 0, we obtain )( 


PB l i ae er, a 
= — En 3 1 
= RA) RE—3E*+E°) (Fig.8) (14), u 
o being measured from the intermediate support B to Ra Ra Re 
the right. Ag 
From the above analyses, it can be seen that Fig. 9 


the method here presented is a sort of “the method 
of initial conditions”, provided that statically-indeterminate quantities are known. The present 
method might be applicable to the analysis of the rahmen and the arch structures. 


However, when the system becomes complicated and its order of indeterminateness is 
very high, the above method would necessitate some modification, and then recourse would be 
given by such asthe SouTHweın’srelaxation method. In fact, simple cases as above may be treated 
by not assuming statically-indeterminate reactions and moments, the method of which is similar 
to that given in the subsequent article. Apart from this line of investigation, we shall proceed 


on the bending of plates. 


3. Bending of Plates 


The first example of the present method applied to plate problems is the case in which a 


rectangular plate is subjected to an eccentric concentrated load, the four sides of which are sup- 


posed to be of simple support. Rectangular coordinatesare taken as shown in Fig. 10, in which both 
sides are denoted by a and b, and the concentrated load P is applied at point A(&, n). 
The deflection w of plate in general satisfies the bi- 
harmonic equation 
E 02 0° \2 0 15 
ae 
in the domain where there are no loads applied. Ifthe plate 
under consideration be divided into two parts (I) and (II) 
by the straight line BC (y = n) parallel to the axis of Oz, 


i 1 \ 
the deflection w, for the first part (I) (Sn <y< n) 7 
takes the form 


w = > [A(x) cho + B&) oshoe + C(x)shoe + D(x) och o] 
n=1 


FEIN 
en. -(16), 
in which for brevity 
nzy nan Fig. 10 
[% ——, , SE 5 
a a 


A(x), B(x), C(x), D(x) being constants depending on the parameter x. Then by the reciprocal 
I Ze 

theorem, the deflection wır for the second part (II) In e4=7 ) is given by interchanging o 

and x in equation (16); that is 

Nn& 


ne Be 


In = 5 [A(o) chx + B(o)xshx + C(o)shx + D(o)x ch x] sin 
n=1 


Equations (16) and (17) satisfy the condition that the sidesx = O and x = aaare of simple support. 
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Now since the function (17) must satisfy the biharmonie equation (15), A(®), - - - , D(x) take 
the form 


Alk) (Ba u chx 
Ba) _ | ds db ba) | abe 
0) GG & sh x 
D(x) dd d d »chx 


where a, a, ... , d, are independent of x. In what follows matrix form will be used exclusively, 
which greatly facilitates our calculations. Then by substitution equations (16) and (17) are 
written in the forms 


u GG % ch x 
S bs bb b& 5 »shx | . nzz 
= sin a FE 
w; ze [ch 0, osh 0, sh o, och o] DIR AUWS en = 
der dad, »chx 
dla 08. 0, ch o 
® b, eh oeshe|.. nzx 90): 
DW = 2 x,25hx,shx,xch x] en EEE RR (20); 
derdede 0) ocho 
and the latter can be rearranged in the form 
rcberd chx | 
> PA er xsha I Nnn® 41 
Wy a [eh 8, esh 0, sh 0, och o] a : ch, j a 
N PERF xchx 
The coefficient matrices consisting of a, a, ..., d, in equations (19) and (21) are of transposed 


matrix each other, and it can be stated that the reciprocity in equations (16) and (17) is the same 
thing as the transposableness of the coefficient matrices. 


The problem to be solved then reduces to determine the coefficient matrices so as to satisfy 
the continuity conditions along the straight line BC where y = n and the supporting conditions 


h b er : : , 
along the sides y = + 7° These conditions will be taken into account successively. 


The continuity conditions are the continuity of deflections, slopes and bending moments 
of the domains (I) and (II) along the straight line y = n, that is oe = x, and as to the shearing 
forces there must be a discontinuity caused by the point load P, which may be represented in the 
form of Fourier sine series. Thus they take in this case the expressions 


u I Bee 
do do ' do? 00° 
’ 2 E =%) 0... RN2R 
uw; OP 2Pa& 1 ELEL ER BELELE. (e ) (22), 
= —- sin —— sin — 
00° 00° n®D u m a a 
n= 


where Dis the flexural rigidity of the plate, i. e. D = E R/12 (1 — »2). The first of equations (22) 
is identically satisfied in virtue of the forms of equations (19) and (20). The second to the fourth 


in equations (22) will be considered. To do this it will be convenient to refer to the successive 
derivatives 


u ch o PAR, Ba eg 4 ch o 
0) 
26 osho s+% 4 a4 +% [28 osho 
® \[,&%&e = 
En ı 62 & el sh o «+2, & 42,6 shoe|' 3), 
[% 
03 


0; och oe eg + 9.64: Bra Ana och o 


j 
Fe 
Ä 
er. 
2 
Br. 
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where e represents one of a, b, c, d. Then after a little rearrangement the continuity conditions 
(22) afford 


a NR a, Q, De 
bb & b, Ay b, bb, 
Rs 2 b, 7 ee (24), 
dd u dd, Or = b; AH 


where for brevity 
KEAE se Aare 
Re 


Itis noted here that, as to the last of equations (22), reference should be made to the identity 
ch’«—sh?x =1. 


The remaining conditions are then those along the sides y= 24 Assuming these sides 
to be of simple support, we have 2 


o:w 
w).-ı=0, I EFT DEE 26), 
wei er (26) 
” 0°Wyr 2 
(od; 6 | a Jule DOREEN DEE (27), 
where A = —_ . The formers of these equations afford respectively 
Babe ch? 0 GT chi 
Drbme, rd, )shA RS AshA x 

— Une 28): 
a a a 5 —-sh} ö Er Dee ca aRc, sh A 2 

„bh 4a d — Ach} Anedseedssad, AchA 


As to the latters of equations (26) and (27), it is preferable, instead of the forms as they stand, 
to write 


2 2 
(n) 8 (er u) et ren): 
0oe=—1 - Jjoe=4 


Referring to the derivatives (23), equations (29) afford 


(d,, d,, dy, d,) = (by, by, b3, b,) eth A, En 
(ar de 64 de) = — (Aa D3, 6, da) cth A 
Thus all the coefficients a,, a, ...., d, can be determined by equations (24), (28) and (30). First, 
substituting equations (30) into equations (24), we have 
I & & 
(1 a Fa Fr a, a th1 @ N 
& 
ey = T th} 0 - 0 
= Se 
\ & & 
N N: Az + 5 2 & zerh 1 
& A 
d, d, d; d, = “T 0 a cth A 0 


Equations (28), consisting of eight equations, must be satisfied simultaneously by a set of certain 


values of @, @dy,. .. ,dy. We see at once that the second, fourth, sixth and eighth in equations (28) 
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are identically satisfied, while the remaining four equations, with the substitution of equations 
(31), afford 


aa,tHi = — (A 2thA—Ath®)), 
3 —gtha=— 7 A(th2—cthA), 

a atha=— TAl—theR), 

+ GthA=—— (2—Athi +4cthA). 


The first and the third of these equations afford the values of coefficients a, and c,, while 
the second and the fourth afford those of a, and c,. Thus we obtain 


ne -{ (A —th2—4th?}), 
9=-—- Kane ee in > (32). 


Bi I (A — cthA—Acth?}) 


It is to be noted that in this case the same value of a, has resulted from the different sets 
of equations, and this would sometimes provide a practical — but important — means for checking 
that the whole system treated is a possible one on the standpoint of mechanics, and that the 
preceding calculations are correct. 


The remaining coefficients q,, a, etc. can at once be determined by substituting equations 
(32) into equations (31). In this way the coefficient matrix in equation (19) has completely been 


1 
determined. Hence the deflection w; for the domain — 2, b<y<nbecomes 
Ip 2 [ch exe ano Eltaretge] 
n= 
—A+thA-+4th24, —thA, —1, 1 ch x 
& —ıu% 0, 1, 0 »shx |. nıaz > 
Er L, —1, A—cthA—Acth?4, cthA shx |" ae 
—1, Dr cth A, 0 xch “| 


and accordingly, by transposing the coefficient matrix, the deflecetion w,; for the domain 


= U< > b becomes 


ws = & [cho,osho,sho,och.o] 
n=1 


—A+th}+Athı, —thA, 'R —1 ch x 
& — th}, 0, — 1, 0 »shx |. nn 
x in . 
4 —]1, IL: A— cthi—4cth?4, cthA sh x ae a > 
1; 0, GEH, 0 x“chx 
where, as before, 
DAR RIEI : DER nzb nam nzy 
ad 4 ag de Duap 9 ; Ale er tea 


We see that the solution obtained in equations (33) and (34) is in accordance with that ob- 
tained by the ordinary method of analysis, and that time and labour is much reduced by the pres- 
ent method. Numerical evaluation of the solution may also be effected efficiently from the 
above forms of solution. 

As a simple extension of the above solution, we take a two-spanned continuous plate sub- 
jected to a concentrated load P. Reference is again made to Fig. 10, in which the axis of Ox is 
supposed to be the supporting line, all the other conditions being assumed to be the same as 
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those of the single-spanned plate above. Let concentrated loads Pa(£) dE be distributed along the 
axis of Ox, gq(&) being a distribution function. The deflection wir in the domain 0 < y< - b due 


to the distributed load is obtained by putting x = Oin equation (34), and by integrating the result 
from O0 to a; that is 


& P®@ xSıl nme 
in = 2 | An 7 dE- [ch 0, 0sh 0, sh 0, och] 
> 0 


x(—A+thi+tb, —thl, —1, 1) sin — N: (36), 


where, for saving space, the parentheses ( ) denote a column vector. 
On the other hand, when a concentrated load P is applied at point A(&, n), the deflection w, 


in.the domain — 9; b<y<nisgiven by equation (33), which for y = 0 takes the value 


P&@ 1 na 
u = ap mi = [eh x, xsh x, sh x, x ch «] 


SR lttnärtälthiAnestun ie Msn Det, et 
a 


Since the deflection along the axis of Ox must vanish, we have 


(wn)y=-o + wDy=-0o—=0; 
from which 


WERE e 


@—thA—Ath?}) | g(&) sin 


sr [ch x, x sh x, sh»,#chx](—A+thA +Ath?A, —th,— 1,1). 


—= sin 


To determine the function g(&) we put 


1d=Dasinı, 
n=1l 


a 
and then 


a 


[ne sin "FE as = cı | u et 


By substitution c, is found to be 


Deere 1 


a oe ep 


in 22h, th — UN)“ 


in virtue of which wır becomes 


BEW. . NT 
wi = DR [cha osho,sh@achol (—A+thA+ AA —thA —1, 1) sin — (9) 
1 


2mD n= 
Thus the wanted deflection w in the domain 0 < y< nis given by the sum 

A RE LT 6») Pe (39). 
Similar consideration will enable us to obtain deflections for the remaining domains — > b<y<0 


andn<y< - db. Furthermore, the three-spanned continuous plate may be treated similarly. 
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The second example is the case of circular plate subjected to an eccentric concentrated load, 
whose outer edge is clamped or simply supported (Fig. 11). The differential equation for deflee- 
tion w, instead of equation (15), is 


ra a ea gt ran 
vr 


y 


and solutions suitable for the present example are 


1 — Ry(, 0) + Rı@%, 0) 608 0 + & Ruf, 0) cos 0 (40), 
where "> 
Ro(&, 0) = Ao(#) + Bu(x) 0? 
+ Co) logo + Du(z) 0*log o, 
R,(&, 0) = Aı(x) oe + Bı(«) 0? 
+ G@) ET" + Dix) elog eo, 
R„(#, 0) = An(#) 0" + Bu(x) 0" *? 
+ Ca) E”" + Due) oT" +? 


x, o representing the dimensionl ss quantities 


EN NE 3 r 


“= — =—., 
Fig. 11 Wi ° a 


(41); 


Let w be denoted by deflection for the domain x < .o< 1, and then deflection w’ for the domain 
0 <o<xisgiven by interchanging o and x in equation (40) in virtue of the reciprocal theorem. 
Hence 


w = Rle%)+ Rex) co8 + DAule, x) wen 2 2 Pe (42). 
Comparing equation (40) with equation (42), we must have 
Au) a4 % A,l&) a @ Anz) u @ 
Bu(x) _ [dA |]1 BA)| _ |b 5|Il» B&)| Ib & a (43 
cal Ta el Ic |" I& ahel [a w| [er Ihre] 
D,%) dı d, D,%) dı ds Da(&) di ds 
where a, 4, ...,d, do not depend on x. It is noted here that terms which become infinite as 0 


tends to zero are excluded, since there is no hole at the center of the plate. Thus equations (40) 
and (42) take the forms 


a 0 u % 
4b, [fi lie 
w = [1, 02, log o, 0? log u Be z ı h 
Il, g 0, o*logo] ae Bi + [e, 0°, o=!, o log o] a 2] cos 9 
dı d, Me 
rs = 
S : ee gr 
N n+2 —n —n+2 1 2 x } 
+2, a a ] Er +] cosnd „ ...(44), 
u % 
4 4 
A SE. „2 Fe hr, 28 
Ww':= rl 02] 1 1 1 1 9, 3 1 1 1 #' 3 
B od, log x + le e] ab od en Gen) Parse i 
x" log x x log x 
un 


" 7 7 n Re 
ab a a u cosn d (45) 


[e.0] 
+ E37 [0" n+2 
lee | Mb Dan 


nn ie En uch u 
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The continuity conditions along the circle e = x are 
(62) . (5) o:w 7710:0% 
O0)  \0 122 Bel. Ge), 
ow ow’ ag el = . 
en) 


in which the right-hand side of the last equation is the expression that the concentrated load P 
is written in terms of Fourier cosine series along the eircle. We shall treat the continuity condi- 
tions (46), and after that boundary conditions along the outer edge (o = 1) will be considered. 


We now have the successive derivatives 


1 aA er &: ELMETFE, 1 0 0 0 
En 0°? 0 28, +& RT TACH 0? 0 1 0 
2 —— 
Be ne I 2, +34 —& 24 logo 7 ee 
3 - 
ee go? log o 0 0 26% 2 | | o*logo ologo logo oT! 
ed eg ') & & e & Fee 1 0750 
e eo at, 3. —& e4 % 0° ei 1 
Dart fer ee = i 
ER [eı € e3 €] o-! 0 6e, Dyes e, 07 0° a 
Er glog o 0 6, —6% —&| | eloge logo ' e° 
1 p® 
1) 
n+2 
00 ° 
„or [eı € es &ı] = 
2 
® —n+2 
d0° und 
e, 6, eg e N 
ne, (n+2)& — ne —(n—2) & 
n(n—1)e, (n+l)(n+9e, n(n +1), ee oe 


nn) (n— Ye nn +l)(a+9), —nn+D(na +9, —n(ın —1)(n—2)e, 


0" D* —1 og —2 o" —B3 
0" +2 oe" +1 0" DZ 1 
x Oo" 0 n—1 0 n—2 0 n—3 a 
ai +2 Da +1 Pag em 
where e represents one of a, b,..., d’, and the products of matrices should be taken to consist 


of diagonal elements alone. 
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Referring to the successive derivatives above, the continuity conditions (46) afford, with 


equations (44) and (45), the following: j 
Re De a ned a a a; 
b b g,—a b bi bs NER > bi bs az b; 
= ' ' = & " = & 47 ’ 
ae 0 ®  & 0 rs ! [41 & 0 —nm+)) ( ) 
di d, ©. 0 ll 20.0 ni ZELL 0 
P« 


where for brevity & = 8zD' 


Equations (47) are the continuity conditions, apart from boundary conditions along the 
outer edge oe=1. Inrespective matrices three constants still remain undetermined, and they will 
become determinate when the outer boundary conditions are taken into consideration. 


We take the outer edge to be of clamped support, which is expressed by the equations 


dw 
(w),-ı=0, (Gr 1 I Be (48). 
Equations (48) give 
y+b=0, 2b, +4 +d4=0; 
++ =0, +3 +d=0;0 2. . (9 
+ +a+K=0, na+n+Y)b nd — nd —0 | 
where n=1,2,... Substituting equations (47) into equations (49), the coefficient matrices 


become determinate. 
Thus the wanted deflections w (for x < oe < 1) and w’ (for 0 < oe < x) become respectively 


1 1 1 1 
a te 
| 1 1 
w= «[1, 0%, log g, o*log e] +a[e 0%, 7", eloge] 1 cos 0 
0 1 #2 0 —— I 
2 
1 0 —2 0 
men 1 
n—]1l n 
1 1 
Se: : 5 n a+l x 
72% 0,872, 9 1 cosn® . . (50), 
Br rer. ER Eu Eh yn+2 
n(n + 1) 
L 0 
n(n — 5 
1 RR; 
EEE 
Bro] 2 2 * —1 1 0 —2 „3 
el] 1 +0 [o, 0°] 1 i cos d 
D% — Do 1 0 log % jr 57 _— 5 0 „1 
x" log x x log x 
- Im 
1 1 1 
-I- 54 n n-+2 un n 0 n(n — 1) RR 
NE 1 1 1 cosnd (51), 
n n.+-l:ysnaT) : es 


. na Ar 


u a. 


AA 


er 
P 
- 
“r 
“ 
we. 
er 
3 
e 
> 

r 
Ei, 
ar 
ri 
“ 
« ’ 
3 

Fr 

+ 
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= 


where as before 


a a 
8nD’ ei, ER 


Are, forms of expression for equations (50) and (51) are not so complicated, and they take 
the forms 


ende, n : 1 4 r 
Gele ES NET) 
£ P« 1 1 
EP + [a9 5, @- E45 ent + 2elogelcos6 
ProiieQ n 1 x n 1 # n 
TEETNDR al ram ke een re 
x en 1 As 
pn elle Ge | 
Ban Pr a® 3 1 
er, ere)lgr + zit) A—R) 
P.a3 1 
4 ei» ze) + ger +2ulogr cos O 
Pe Si 1 On 1 0° n+2 
8=D een ee 
a Frl Sehe! —n+2 
aa TE BR | os 0 0(<eo<x). 


We see that the above result is the same as that obtained by the ordinary procedure [5]. 
If, instead of equations (48), we take 


ww 10w . 1 dw 
FAT Er | al Tr @ 


we then have the case in which the outer edge of the circular plate is simply supported; » being 
Poısson’s ratio of the plate. The rest of this case is here omitted. 


In conclusion, it is added that the present report is confined to the cases whose solutions are 
known otherwise. But several additional examples, which are more complicated and seem their 
solutions not to have been given owing to the complexity, have been worked out, and they will 
appear some other day. 
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Eine geometrisch nichtlineare Schalentheorie 


Von D. RüpIGER 


Unter Verwendung allgemeiner Koordinaten wird eine geometrisch nichtlineare Theorie elastischer 
Schalen aufgestellt. Aus ihren Grundgleichungen kann man weiterhin die linearen Differentialgleichungen 
der Stabilitäts- und der Spannungstheorie ableiten. 


In eurvilinar coordinates, a geometrically nonlinear theory of elastic shells is derived, and the basic 
equations are given. From them we obtain the linear differential equations of the theory of stability and the 
theory of infinitesimal displacements. 


C TIOMOINBIO 061MHX KOOPAHHAT BbIBOJUTCH TEOMETPHUYECKH HelIuHeäHaA TeopHA YNACTHYHBIX 
00010Yek U YKA3bIBAIOTCH UX OCHOBHLIE YPaBHeHHA, U3 KOTOPBIX MOFKHO IIOJIYYHTb JIHHeÄHBLIE 
andepepenumasbHpie yPaBHeHHuA TEOPuUm YCTOüyHBOCcTH U Hanpsukennä. 


1. Einleitung 


Unter den Voraussetzungen der BErnouLuischen Hypothese und der Unveränderlichkeit 
der Schalendicke h bei der Deformation wird im Sinne der Theorie flacher Schalen großer Form- 
änderungen [1] eine geometrisch nichtlineare Schalentheorie aufgestellt. Durch die Anwendung 
schiefwinkliger Koordinaten gelingt es, die Grundgleichungen in allgemeiner Form abzuleiten!). 


Zur geometrischen Beschreibung der Schalenmittelfläche und des Raumes innerhalb der 
oberen und unteren Schalenleibung werden einige Ergebnisse aus der Theorie allgemeiner Koordi- 
naten und der Flächentheorie verwendet. Die rechtwinklig kartesischen Koordinaten eines Raum- 


punktes sind x bzw. x*, xP, x, wobei die lateinischen Buchstaben an Stelle der Zahlen 1.25 
stehen. Die drei krummlinigen Koordinaten eines Raumpunktes werden mit x* bzw. x#, x”, 2 
bezeichnet. Die griechischen Buchstaben }, u, v, o kennzejchnen die Zahlen 1, 2, 3. Die karte- 
sischen Koordinaten x‘ sind differenzierbare Funktionen der krummlinigen Koordinaten x” und 


umgekehrt. Durch partielle Differentiation 


U) %) 
a On; zu On 


erhält man die Tangenten- und Gradientenvektoren 
ch = a s 2 = x k EEE, ae (1,)- 
Für ihre partiellen Ableitungen wird 
k 
CA = Au s Chu = Cu lin ran ae ee ee (1,) 


geschrieben. Die Übertragung des KRoNEckersymbols 6;, = ö'* in die allgemeinen Koordinaten 
liefert den Maßtensor in kovarianter und kontravarianter Darstellung 


= nk A NESAT 
Nu dr eh Mae eh ra re 


wobei vereinbart wird, daß über gleiche, in einem Produkt oben und unten auftretende Indizes 
zu summieren ist. Die Unterdeterminante Yıı Jda— N? der kovarianten Koordinaten des Maß- 
tensors wil d mit g abgekürzt. Die kovarianten Ableitungen im Raum werden dureh einen langen 
Doppelstrich 

AA 1 : 

A Ins B*u||,, Bill, ee TE EU BIBTT u u, Ba (3,), 


Überschiebungen der kovarianten Ableitungen mit dem Maßtensor durch 


ae 
Ar A|ugPt. 2 58 An ee, A 
gekennzeichnet, Die Curistorrevschen Dreiindizessymbole zweiter Art sind durch die Glei- 
chungen 
a rt ERROR. 
Dr 0 = u.we54 (den Tr Ga (4) 


gegeben. 


) Hinsichtlich der gewählten Bezeichnungsweise vgl. [2]. 


un u BEP UUETEN 


TER mir Fr m TEEN TER FR TAIREN 
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Zur geometrischen Beschreibung der Punkte innerhalb der Schalenleibung dient ein allge- 
meines Koordinatensystem, dessen Koordinate x? = z mit der Richtung des Normalenvektors 
N* = N, der Schalenmittelfläche in jedem Punkte identisch ist. Für dieses System verschwinden 
wegen ON, =&N, =0, dA N*=&N,=0 die gemischten Koordinaten des Maßtensors im 


Raum = 9° = 5-00, und es ist 9? = Go, le 


2. Geometrie der undeformierten Schalenmittelfläche 


Die kartesischen Koordinaten der Schalenmittelfläche vor der Deformation seien x! bzw. 
x*, &P, &2, ihre Flächenkoordinaten x* bzw. x, x”, x, x°, x". Die griechischen Buchstaben «, By; 
ö, 0, t sind im Gegensatz zu ), u, v, o Repräsentanten der Zahlen 1 und 2. Die Tangenten- und 
Gradientenvektoren des Koordinatennetzes der Schalenmittelfläche betragen 


En 0 FÜR: 


Ca = an Ge a ee N. 
Für ihre partiellen Ableitungen wird 
le, var u here a a ee) 
eingeführt. Die kontra- und kovarianten Koordinaten des Maßtensors sind durch 
BR 6 0 VE BR N a ki) 


gegeben. Wird die Determinante der kovarianten Koordinaten des Maßtensors Det a,,; = 
Aıı dgga — (ya Az, mit a bezeichnet, so erhält man für den Normalenvektor' 


N, = — &pg ® € LE or ee a u ar 9 


ya 


Darin ist &,, das bekannte Permutationssymbol. Die Koeffizienten der zweiten und 
dritten Fundamentalform lauten 


Da, Web, ee! 8) 


Ci Day; cd, bi, ern, x 


& 


Für die Determinante Det b,; = bj1 dag — dia dzı wird die Abkürzung b verwendet. Die ko- 
und kontravarianten Differentiationen auf der undeformierten Schalenmittelfläche werden durch 
einen langen senkrechten Strich 


A* 


B> Bei, : AsI|# ae A® 


OT ee eek) 
gekennzeichnet. Zwischen zweiten kovarianten Ableitungen besteht die Beziehung 
Ed U RE U) 2 Ve NH 
Die CHRisTtorreLsymbole zweiter Art der Schalenmittelfläche lauten 


o \ 1 
SR EN &Ö 
5% u a Ch Chy = 5) (Bi (Ay, - Apgsıy > Apyj6) re: (LUD): 


3. Geometrie des undeformierten Schalenraumes 


Vor der Deformation kann der Ortsvektor x" der Schalenpunkte im Raum innerhalb der 
Schalenleibungen in der Form 
A ed > RE En 
dargestellt werden. 
Daraus folgen durch Differentiation unter Beachtung von N* 
Schalentensor 


ok 


I: B Zi 5 
«= -5% mit dem 


NER ET (18) 


die Tangentenvektoren zu 
a RR ER AA). 


a a a on die G! 
Damit können die für eine weitere ee Rech benötigten Koordin \ 
DeBISTorsRLEyEnNIolE zweiter Art ermittelt werden. Für diese gilt (2), Bl 4 
Gap = Asp — 22 bag + Ptap, Da aeg 
=ah+z2ur+@deh, ‘0... 2 I 
& N WR ist ri or .% 

gl, Ia=f?=0) m F 
| her = 

pr = Per — 205, BE Bl Bin: 

T!yp = bag — 2 Cap » 


u | 


n r 
u EERTER HOSE FEER, m: 


Tas = —-b —25 — 2% b2, 

Pa=Ta=Ta=0 - u u = 
Die kontra- und kovarianten Ableitungen des Vin ivertnne im Schalenraum ergeben : sich 
dann zu E 


I ar 
a De 5 


ve] = ve — be# VE + z(—3cP VE + 208 Ver — BEP vr) + 
| +2(—6berch VE + 3.07 Vel, — cf], VO — BPrbi]|, V?), 
v=VP—%Vv+z@8V,’ +2 v—&V)+ 
+2(8c$ V,P +75), V?+c,V2— 03 V3), Y 
v3||® = V5|@ 4 58 vr + z(c® Vr + 252 Vep) + 22c2 (8 VE+ 3 VEp), an). # 
volle = V|e +58 vr + 2802 vr + 25% Vep) + 230% (255 V2 + Vep), | 
vh=V).—K Paz Pr -+ EA ENTE 
Ve|ß + Volle = Velß + Vele — 2 beB VL 22(b% Ver + 58 Ver — 5 Vr— 3cP V3) + 
+ 223 (c$ VPlr + .c# Velr — cP|,Vr —4c$ bPr v9) 


. 


Mit ihnen kann man die Ausdrücke 
1 
> „IP vrle + V,]P Velo 
= Vrev,— (54 vr? +55 vr) VB +? VE VS + z[2(b Vr|ß + 5 Verl) Vo — 
— 2(c2 VrP + ch Vr%) VE + Ab&cPr ve vs —2(bE 5 + b£ 5°) vr? Ve] + 
+ 2 [3 (e$ Vr|P +c$ Ve) V„@+ 45552 Vrjr Vo + 108 cdr vs ve — Dede v, vo + 


1 1 
tz @ var vNBV 2 ar aD VW 


1 
— 5 68j° VrP +55 VD BE V.— Zd8 (ed VrP + Vrp) Vo — 


—3 (4 b5 + ch b2) vr[® v3 —4 ae bi Ach) vr ie Be EEE (18,) 
und 


1 
> (V3]l* v5]? + v2] V3]%) 
= Vve|* ve]? + (by v2]? +5 Ve|%) Vr + 5558 Vr VO + z[2(b% VS|P + 52 vsje) var + 
+2 = b5 + 65.b2) V? V2lr +2ce$ VB]? + ch Ve|®) vr + 20% 08 + 52ch) Vr Veh + 
+21; um? v2? +55 ve|$) cd vr +4(bich +52) VE Ver + 3ckch vr Ve + 
+3 E2 cs + by. 63) v7 V2|P +3 (5 V2]P +c$ Vope) var + 
7 
+5 b6 + BE b5)cx VO Ve HABEHB vewvalol. . . oo. (18,) 


bilden, die für die weitere Rechnung benötigt werden. 
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4. Geometrie der deformierten Schalenmittelfläche 


Bezeichnen V* die kartesischen Koordinaten des Vergeiebungsvehtgrs der Schalenmittel- 
' fläche, dann ist deren Ortsvektor #* durch 


EEE 257 un (9) 


gegeben?). Unter Beachtung der Transformationsregel 


Raw ınw 


folgt a | | 
Pen Lo VvieNnsWwe,.. m (20). 
Durch Differentiation erhält man daraus die Tangentenvektoren 
a EN NER VE: al 


die mit den Abkürzungen 


o 1 C) ° o - 1 5 er A 
As Ei Do (Ve + V?) — be? W A Aus == oJ (Vals ir Vale) == by Ww, F (22) 


As va, HE W 


die kovarianten Koordinaten 


Gy = Aug + 2Äse + (Va — 5, W) (Ve —de, W) + (bay V’+ Wie) (bes VE + Wi) (23) 


des Maßtensors der deformierten Schalenmittelfläche und die Determinante 


2; 


o 


Deta,g=a=a[l1 +2A: +2Ä, AH — 24544 4 
+ -EWIE EWR ERICH 


festlegen. Die kontravarianten Koordinaten des Maßtensors, die Gradientenvektoren und der 
Normalenvektor der deformierten Schalenmittelfläche ergeben sich zu 


GB — a _2 Asp 1 4 As Avp — (Vol — dee W) (| bEW) — 
— (b2 Ve + WI) (6 Ve + WI) +: -, 
2-8 +8[-(W,—b£ W) +2A3 (v,— 55 W) — 
— (1,35 W) (5 W) — li V7 + WI) EV +WI+HL, . os). 
+ N# [56 v7 + wP— 2A (e vr +W@]+---, 


a Del Bdız D ken a N A ala Hl al a a ia a an u 0, 


a A 


—— 1 o o o © & 

N, = N [15 0, 9° + Wi (0 9° + we] — 

— [day + WE + WP) (Vo bap WI + 
| Für den Haupttensor gilt 


| Dep — das + Was Ca W + 5 Vo + 25 Vla+ 5, + bat tat» | 
EZ ı wp? +3? w— 8 Ver — 8 Ver HP VW + Hr, 026). 
b NN och ae Kl 
b5 = WE +EW— 8 VL, +5 V, +, V 4. =bhtn,d4e 


Damit werden die Determinante 
Detb,, =b = 5 [53 (w wie—cı W + 52), Ve + 202 V’,)— 52 WI + 
+. W—b2b, V—2 WL]+-, 
b b b \70 au ö (yirle ‚our b° vr J.h° v' je: 
ed kV 53 W):.+ b5 (we —3W+ = 2 T \% 
a a 


2 ” b 
pe WIE +. W—2b,, ”’ —2c vl a a v2) 


>07 


2) Vgl. hierzu [6]. 


14 


Ze „ BG ine h Kaas Fee 
i ZN as 3 


een 1 5.05 VPr—b 
= ange er 


= ELLE A TEE R ER ale 
eb Aa) 


ermittelt. Dabei stehen die. 


bungen linearen Glieder. Die kovarianten Ah ı 
_ auf der deformierten Schalenmittelfläche cas 


ud ur RTumätg+ Puär, £ 
BR, = tn, u 


die CHRISTOFFELsymbole zweiter Art für die deformierte Schalenmittelfläche. Sieh aben die 


Werte T | 
c= © © EIER 1 © - a N E > « 
| Dias Tao + da (5, "+ WI) + (Vlat FW > 
7 2a 
i VW) + | 
Mit y o ° ao o Bi 
, Vlaa= Vlag + baobz V’ — bas b; v’ \ 
erhält man etwas kürzer : a 


Prag Peap + bag (65 7’ + WI) + Vlas—d5 82 0, 5 WEW + DEWID + --- (2). 


5. Geometrie des deformierten Schalenraumes 


Wenn man die Koordinate in Normalenrichtung der deformierten Schalenmittelfläche mit 
z bezeichnet, gilt für den Ortsvektor jedes Punktes innerhalb des deformierten Schalenraumes 


SPAN ... 0.0 EEE (33). 


Nimmt man an, daß für die endlichen Deformationen z = z gesetzt werden kann, so ergeben 
sich mit dem Schalentensor 


Fl = 8° — z° 


der deformierten Schale die Tangentenvektoren zu 


Ef nn Me Eee (35). 


Aus diesen können die ko- und kontravarianten Koordinaten des Maßtensors im deformier- 
ten Schalenraum 


IT lg — 22 bg + 2iyp, | 


gyP=aP+zIbB +23 CHR 
sowie 


YDet 5; = 9 = Yalı-z3 +22) E (37) 


gebildet werden. 


D. Rüpıser, Eine geometrisch nichtlineare Schalentheorie 203 


6. Deformationstensor für endliche Verschiebungen 


Mit den kontravarianten Koordinaten V* des Verschiebungsvektors lauten die kontra- 
varianten Koordinaten des Deformationstensors für endliche Verschiebungen bekanntlich 


1 1 
D** — DE VA Hr yore = (v,|? vP + vr vr) NL ERMEREE, (38). 


Die in der Theorie der annähernd zweidimensionalen Kontinua benötigten Flächenkoor- 
dinaten betragen danach 


& 1 & & 1 & & 1 | 
Dr = |vel + ve + ev + vn + le voll + Valle ven] @9. 


7. Elastizitätsgesetz 


Es wird angenommen, daß zwischen dem Deformationstensor (39) und dem Spannungs- 
tensor der deformierten Schale das Hooxzsche Gesetz besteht. Dieses hat für ein homogenes 
isotropes Medium mit dem Schubelastizitätsmodul G und der Querdehnungszahl » die Form 


#26 (D* u 7 ge# Dy) a Re (40,) 


und geht mit der zweiten Gleichung (15) in 


y 


8 _9G DP+ — a’ Dy+- (2bPz2+3c#)D| .... (40) 


1-—» 
über. 


8. Bernoullische Hypothese 
Die Verschiebungen V’ jedes Schalenpunktes sollen so sein, daß mit 
VEN T VO DNM. . ee (41) 


alle Punkte, die sich vor der Deformation in einer Normalen zur undeformierten Schalenmittel- 
fläche befinden, auch nach der Deformation in einer Geraden senkrecht zur deformierten Schalen- 
mittelfläche liegen. Bei Beachtung der Transformationsregeln 


v_el, won 
erhält man dann für die Verschiebungen V*, V? des Schalenraumes 


vV’=-eV= v* + z [2 w[° 2. (v? %__ p%B w) % r 3 w|,)] % 
—+ z2 b}; IE W ß a (ver ig pP? w) (b? v, Pr w|,)] 


le, 
V-N,Vv=-W-—z7 (9 V,+ Wl) (d$ V’ + WI?) 


Die Gleichungen (42) ermöglichen es, die Flächenkoordinaten (39) des Deformationstensors 


als Funktion des Abstandes von der Schalenmittelfläche in der Form 


pre A ı Berrz(Aerı br) +zlAeı Be), | 


< s “= EEE | 1.2 (43) 
D} = D*Eg,5 = As + Ei 12 (As 7 B2) ur 2 (A: % B:) 
anzugeben, wobei die in (43) verwendeten Abkürzungen folgende Bedeutung haben: 
Aa® — _ yes 4 pe Vor + 98 vor — BEP Pr — 3? W, 
1 o 

ME A 25, AR, | 

2 3 ra B a 3 m 2 .B1y « B Ta|y &% „By Sy 3 “B vr (44,), 
ABP = zw’ BuW N HE (c& ver + ch ver) —6 bc?’ W ai? | 


2 2 10 t be 
A? = Ayß ArELS Dyg Ar “ + CyB A 
14* 


Edz cry de 
rt We | | 


Ber — 2 wis Wi? en ww + ir ,(& ee + eu ) 


2 


2 
B; =” Axp BP — 2 bp Br# + CB BP 


es 2 T on j ea. a Er u ". wea = 0 


Ne LIKE Br 


Bl an +z° a 


ek: ee wir Hager 1 (we mw) vr 
4 Valzaaa TE w er va ee a 
a shi P 


© wu A 


us WE Lam vr Hp a ER 


Ic E- 
2% (dt W|® + 58 WI B Vrsp EEE ie (3 ve +» Ba ! er a 


3 2 I B er 
re BB + ch »» vor W—3 (BR +0 8%) vo Ww= ; ji 
5 Ey 
Be, PW— 2% os pr>], vr Wo +) VW + 7 

IA 

3 o © o o o o, o 
a ne -- 63. (05. V°|® + D8 veje) Vo? + Be 
3 o © o o © . 0) | ° ® 
+5 la ver res VD, 420508 Ve? Depp + Dar (05 Vo]? +28 Velo) V.+ Ri 
‚| 9 1 0) © © R n 2 
+56, be + bEb;| + Va ehn Vz osjrot], Klo 
1 © o i 
el I Pr + erg, ver E 
a 
una. ke 177% ; 
+75 +55) VV aaa? | 
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9. Schnittgrößen als Funktion der Verschiebung der Schalenmittelfläche 


Die in der deformierten Schalenmittelfläche wirksamen Längskräfte und Momente werden 
analog zur linearen Theorie ([7]) durch die Gleichungen 


2 _ na 
no = | Vorne, ur | Wr nee NT. (45) 
—hj2 —hl2 


definiert. Die kontravarianten Koordinaten des Spannungstensors kann man mit Hilfe des 
Elastizitätsgesetzes (40,) und des Deformationstensors (43) als Funktion des Abstandes von der 
deformierten Schalenmittelfläche darstellen: 


y 
1» 


rl Är tier a? (Ay + BY) + 
+e[äAr + Ber 4 t larläy 4 Br) + 2000 (ä; + Byj)| 


+z|Ae 4 Ber le (Az + Br) + 208° (Ar + Br) +3 ® (Ar + 3) (46). 


Bei Beachtung von (26), (27), (28) und (29) ergeben sich somit für die Schnittgrößen (45) die 
Beziehungen 


y 
1—»v 


NaB_ Tan IR VORDERE IE 7 30 B aß Pv 
N =26nlA Heat Br 


„h2 
Sn 


o ; b = 
s ige Re Arr (by 1% 2 +47 8%, +@,h)+ 


2 1 o b 1 
a ei ich 
5 b 
[aa 4 Ar 0 — at g) 4A} (A Une Bis Ve 
ö b 
ge BY + BY (b*® — a®? d$) —. By (er —or nee 


N (rt rar) + Ä} (rer par 2 0er a + .: | (47,), 


Br Gh® 7 8° 1 o f o 
Mr ER (de Aero Ha + Br Aa + 


A u Hy (na ü 
ai —— (@’ 4; + Är (br? — a®P H8) + a®? BY + Bye — aa? 93) — 


a N 


10. Gleichgewiehtsbedingungen an der deformierten Sehalenmittelfläche 


Die fünf unabhängigen Gleichgewichtsbedingungen eines Schalenelementes können aus 
der linearen Theorie?) übernommen werden, wenn man beachtet, daß die Schnittgrößen N®B, 


Q* und M*® in der deformierten Schalenmittelfläche wirken, und daher auch die kovarianten 
Ableitungen auf der deformierten Schalenmittelfläche zu bilden sind. Man erhält mit den 


CHRISTOFFELSymbolen (32,) aus 


Ned rt, MT,— 0 =0..(8) 


s) Siehe [7], GIn. (10.4.6) bis (10.4.1), S. 383. 
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bei Beschränkung auf Glieder bis zur zweiten Ordnung die fünf Gleichgewichtsbedingungen 


N, — 60° a (Ve, — bx w,— bil, Ww) N»? + 
+ [VPla, + (bay 9 — ban BB) V? + du, WI — 5% Wi. — be WI, — 52], WIN’ — 


— (WE + 8 W — 53 VP], + 58V” |. + bil, 7) Q° + pP -0 (4) 
Q* 5 bys N ar (VB: Ww,— DAL, W) Q’ + 
+ [Was — cs W + Bio Vy + dk Vola + b5 Voll N“? + p® =0_(49,), 


Mer, —QP + (Ve, — 52 W,— 55, W) MP + 
+ V Play + (Bay 5 — das 59) V’+ ba, WI — by Wie — 55 WI, — bu, WIM” —=0 (499) 


der deformierten Schale. Die Belastungskomponenten p, p® stehen dabei mit den Belastungs- 


komponenten p?, p? der undeformierten Schale in Zusammenhang. So gilt zum Beispiel für eine 


hydrostatische Belastung pP =p = 0, P =p?=4q. 

Die Gleichungen (47) und die Gleichgewichtsbedingungen (49) bilden ein nichtlineares 
partielles Differentialgleichungssystem für die dreizehn unbekannten Funktionen N*?, M*?, Q*, 
v, W. 

11. Quasiniehtlineare Theorie 

Ebenso wie in der Theorie der elastischen Stäbe kann man auch für die zweidimensionalen 

elastischen Kontinua ein lineares Differentialgleichungssystem für die unbekannten Funktionen 


herleiten. Bezeichnet man die Vorspanngrößen, die für die undeformierte Schale ermittelt 
wurden, mit 


E= 
D*, BD, N*®, Mm, Ö? ’ 
so genügen diese den Gleichgewichtsbedingungen 

+ + + + + 

NP — 0 +0, Ola + bs N +B—0, Mr, —0=0.. 60). 


Bei dem Übergang in die durch Vv*, W gekennzeichnete Lage der Schalenmittelfläche möge die 
Asa 2: 
Belastung p”, p? die Werte 


P=-p#+P, P=eprP,.. mL Deere 


annehmen. Gleichzeitig werden die Vorspanngrößen verzerrt und gehen in 


ee 
über. Die Schnittgrößen der deformierten Schale bestehen dann aus den allein durch die De- 
formation bedingten Anteilen N° ?, M*®, Q* und aus den mit den Vorspannungen zusammen- 


3 ; T + En 
hängenden Anteilen N*P, MP, 0%, so daß 


7 s + 
Nt=NP4ıNe, mMme-MmeıMme, W-ıie ..,. 


. e # 3 e + + + + 

gilt. Zur Festlegung des Zusammenhanges zwischen den Größen N*®, M*P, Q* und N®#, Me? 0° 

2 . 5 x BE # rm . . * .. . 

ist im Rahmen der quasi-nichtlinearen Theorie eine Hypothese einzuführen. Nimmt man an 
Sp x x Se 2 hr MER . R 

daß die Produkte aus den Vorspanngrößen und dem Flächenelement bei der Verschiebung kon- 


stant bleiben, so erhält man aus NEL d= Ne® df, Mm dj= a? df und 0% d= Ö° df mit 
df = Ya da! dx® und df = Ya da! da? die gesuchten Beziehungen ([8]) 
be VE her = anien, 


a 


T ln RE: 
Msn — nie = (1 Ay) air, EN ee 


T 


= 76 -(-aNb 


ir 
jr2 
be 
Di 
“ 
\r 


SE TA 
IR 


TUNER ETERTETDN 


TIEREN EN 
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‘ Die Gleichgewichtsbedingungen (49) nehmen dann mit (52), (53), (51) und bei Vernach- 
lässigung der quadratischen Glieder die Form 


N® A — 580% + [VPlay + (day 55 — Bus BE) VP + Buy WI? — 58 We — BE WI, — 
— 52, W] Ne? — [Wie + & Wr Ve], + 08. Vj, + 02, vrlöe + pP + pP As =0, 
Q* x Zw Daß N a (Wias — Cap w 2; bis % +b V,l on b} v,LNe? en, D° 43 u 
MP —QP + [VPlay + OayB5—ba352) V? + bu, WI — 58 W|,—b2 W|,— 02, W] 7° =0 
an. Zusammen mit einem linearen Elastizitätsgesetz?) 
NeP—ne(V,VeW, Mer me VE Var WR (55) 


liegen dreizehn lineare partielle Differentialgleichungen vor. Sie sind die Grundgleichungen der 
bekannten linearen Stabilitätstheorie. 


(54) 


12. Lineare Theorie 


Durch Vernachlässigung sämtlicher nichtlinearer Ausdrücke werden aus (47) und (49) 
die Grundgleichungen der linearen Schalentheorie gewonnen. Sie lauten unter Verzicht auf den 
Querstrich 

o v o h2 2 1 
N®P=2Gh 1A + Fa a Ar + een — A? (bi + ab) + 


o b y 2 2 Fi 
7 An 5 + )+ rang (’ Ar — a) + 


s ....(56). 
Gh3 


rn y 1 o y o 
Mr= [A + Ar + det) 


N, —50 + pi 0; Oo N MA 0 0). 


In diesen Gleichungen sind jedoch noch Terme enthalten, die dieselbe Größenordnung wie die 
von vornherein vernachlässigten Querkraftdeformationen haben und somit zu streichen sind [9]. 
Es ergeben sich dann aus (56) die Beziehungen 


Be A&B vo oa, FT 3 are __ 1 pa uar ywieß__ 
N =26nÄ + er i4,-55W 5 WI +, W| 
Yy o © & 
le WW, (58), 
3 /o o 
mer (we? ar Wir 
6 —y F 


die zusammen mit den Gleichgewichtsbedingungen (57) ein lineares partielles Differential- 
gleichungssystem für die zwölf unbekannten Funktionen Nager My Ve una WW 
bilden. 
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Eine Begründung für die Anwendung der Laplace-Transformation 
auf nichtlineare Differentialgleichungen 


Von H. LEIPHOLZ 


. Für ein abgeschlossenes Intervall der unabhängigen Veränderlichen wird die Lösung von nichtlinearen 
Differentialgleichungen mittels Laplace-Transformation unter bestimmten Voraussetzungen begründet. 
Durch Zurückführung auf eine gleichwertige Integralgleichung wird die Konvergenz des verwendeten Itera- 
tionsverfahrens gezeigt und es wird eine Fehlerabschätzung gegeben. 

The paper deals with the solution of non-linear differential-equations in the interval [0, t] of the inde- 
pendent variable t, by the application of operational methods. Under well defined hypothesis convergence 
of successive approximation, uniqueness of solution amd error estimation ‚results from the consideration of 
an integral equation which is equivalent to the original differential-equation. 

Ina cayyası 3aMKHYTOTO UHTEPBAIA He3aBuCHMEIX TIePeMEHHBIX IIPM OTPEeleJIöHHBIX YCJIO- 
BUAX NAöTCAH pemenne HemHelHLIX AuhhepenmmanbHbIXx YPaBHeHnÜ IIOCPeACTBOM IIpeoöpa30- 
BaHnun Jlamnaca. ÜBeneHnneM K BKBUBAJIEHTHOMY HHTETPAJIBHOMy YPaBHeHNIo TIOKA3bIBAeTcH 
CXOANMOCTL HCIOJIB30OBAHHOTO HTEPAIMOHHOTO METONA; YKAaabIBACTCH ONEHKA HOTPEIIHOCTH. 


1. Einleitung 


Neben anderen haben die Fachgebiete Schwingungs- und Regelungstechnik das Interesse 
an der Lösung von nichtlinearen Differentialgleichungen stark in den Vordergrund gerückt. Dabei 
wurde von verschiedenen Seiten der Vorschlag gemacht, auch bei nichtlinearen Problemen für 
‚das Aufsuchen von Lösungen die LarLAce-Transformation zu verwenden. So liegt von L. A. 
Pıpzs [1] ein umfassender Bericht über ein solches Vorgehen vor, und in letzter Zeit hat sich 
P. J. NowAcki [2] mit der iterativen Lösung eines Anfangswertproblemes unter Zuhilfenahme 
der LArrAcH-Transformation beschäftigt. 

In allen Fällen ist aber die Frage nach der Zulässigkeit der LaPpLacr-Transformation, der 
Konvergenz des im Bildbereich sich darstellenden Iterationsverfahrens nach einer eindeutigen 
Lösung hin und der Fehlerabschätzung für eine iterativ erhaltene Näherungslösung offen geblieben. 

Es soll versucht werden, im folgenden auf diese Fragen eine Antwort zu geben, womit die 
Anwendung der LAPLAcE-Transformation erst ihre rechte Begründung fände. 

Dabei soll für die weiteren Überlegungen die von P. J. Nowackt ausgesprochene Problem- 
stellung zugrunde gelegt werden, da sich an ihr alles Wesentliche zeigen läßt. 


2. Problemstellung 


Gegeben sei die nichtlineare Differentialgleichung 


De) = fa, DIE)  »#, > arm Teak 
mit den Anfangsbedingungen 
20) = ine) vo a ar Brube 
Dabei ist D ein linearer Differentialoperator n-ter Ordnung, (n > 2), mit konstanten Koeffizien- 
ten. Es ist /(&, x) eine nichtlineare Funktion ihrer Argumente, die keine linearen Glieder enthält, 
und es ist z(l) eine Störfunktion. 
Unter Verwendung der von W. W. SoLODOWNIKOW [3] benutzten Schreibweise erhält man 


durch formale Anwendung der LarLAacr-Transformation und Einführung einer Iterationsvor- 
schrift die Beziehung: 


Dis) Kr) = Mal) FR U EI HZ) (3), 
woraus sich für die (k ++ 1)-te Näherung des Iterationsverfahrens 


ER ı SMa:6) + {fe 2} + Z6S) 4 
ergibt. 

M „.(s) ist dabei ein Polynom in s, das nach Anwendung der LArLacz-Transformation aus 
den vorgeschriebenen Anfangsbedingungen der ursprünglichen Differentialgleichung folst. 


IE DL 4 a Te 1 DINTESEET RETTET 
DEN 


FR 
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3. Voraussetzungen 
a) Es sei z(l) eine im Intervall O<t<T stetige, LAPLACH-transformierbare und damit 
beschränkte Funktion. 

...,bJEs sei f(x, &,) eine im Intervall A< x,%,< B stetige (und damit beschränkte Funk- 
tion ihrer Argumente x, und £;. 

Dabei sind A und B passend gewählte Schranken, die ihre Berechtigung darin finden, daß 
sich alle x, und x, als stetige und damit beschränkte Funktionen im abgeschlossenen Intervall 
0=S1=T ergeben, was in Ziffer 4 gezeigt wird. 

c) Es genüge /(x,, &,) in dem in Frage kommenden Intervall einer Lirscutrzbedingung: 


Me; &) — fı S)]| S Ni Ba HN. 2 2,0, 


wobei N, und N, passend gewählte Konstanten sind. 
d) Außerdem genüge finO <t< T der Bedingung 


IR 7er as a le a 
für beliebige Konstanten K,, K,, &, ß und für passend gewählte Konstanten K, und o («, ß und o 
positiv). 
4. Nachweis für Anwendbarkeit der Laplace-Transformation 


Eine Funktion x(?) ist sicher einseitig LAPLAcH-transformierbar, wenn gilt: 


I a An). 
«| <Ket für t>0; co nicht negativ 
Es soll gezeigt werden, daß allex,(r=(0,...,k-+ 1), die sich aus der Iterationsvorschrift 
DER er %) +2), | ei 
N (0er el. na) 


ergeben, dieser Bedingung (7) genügen, wenn man nur annimmt, daß man die x, aus dem metri- 
schen, vollständigen Raum Rt der nO<st<T stetigen, stetig differenzierbaren Funktionen 


auswählt. 

Beweis: FürinO sts T stetiges und stetig differenzierbares x, ist die rechte Seite von (8) 
eine stetige Funktion von t; denn z(f) ist dies nach Voraussetzung und f(z;(l), &,()) ist dies eben- 
falls als stetige Funktion der stetigen Funktionen x, und &,. Für fest gewähltes x, hat man 
daher an Stelle von (8): 


Dax +1) = YO, h. ans 5b Butler. SR 


A rel ....n— 1) 


Das ist eine inhomogene, lineare Differentialgleichung mit konstanten Koeffizienten für die 
Funktion 2; ;ı() mit n O<sSt<sT stetiger Störfunktion y(f) und Anfangsbedingungen c.. 
Nach der Theorie dieser Differentialgleichungen [4] erhält man für x. +1(f) eine im Intervall 
0 <t<T eindeutige, stetige, n-mal stetig differenzierbare Lösung. Durch Schluß von k+1 
auf k + 2 folgt aber, daß die Iterationsvorschrift (8) nicht aus dem Raum Rt herausführt und 
daher in diesem unbeschränkt ausführbar ist. 

Wenn aber jedes x„(r =0,...,k + 1;k + 1— oo) dem Raum Rt angehört, so sind sowohl 
x, als auch &, im betrachteten {-Intervall beschränkt und genügen dort den Bedingungen: 


Ben heuer, 

EIS ER, MM, G=V...,‚k+137k+18)%, PB positiv 
Folglich ist das vorher benutzte B = Max (K, e*T, K,e??) und ebenso läßt sich aus der Be- 
schränktheit von x, und x, die untere, noch unbestimmt gelassene Schranke A festlegen. 


Es ist aber die Lösung von (9) sogar n-mal stetig differenzierbar. Also gilttinOs!sT 
für passend gewähltes y, auch noch: 


We‘ 


Kae, lU—2, 4,7 0),,, 3%, POSIEIW cd such), 


Da außerhalb des Z-Intervalles über x, und seine Ableitungen nichts vorausgesetzt wird, 
kann man sich alle diese Funktionen dort beliebig festgesetzt denken, insbesondere so, daß die 
Bedingung (7) für jede einzelne dieser Funktionen gerade erfüllt ist. Dann kann man aber auf 
D(&,) für jedes beliebiger = 0,...,k+11; k + 1— 00, die LArLAc#H-Transformation anwenden 
— wegen der Gültigkeit von (7) —, bei vorgeschriebenen Anfangswerten c;. 
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Es kann aber auch auf z(f) die Larzacr-Transformation angewendet werden, nach Vor- 
aussetzung; und es kann auf /(x,(l), &,()) ebenso die Larrace-Transformation angewendet 
werden, wegen (10) und (6) im betrachteten {-Intervall und bei entsprechend passender Fort- 
setzung von x, und &, außerhalb des Intervalles. 

Damit ist aber gezeigt, daß alle Iterationsstufen von (8) unbeschränkt LAPLACE-trans- 
formierbar sind, so daß (3) gültig ist. 


5. Nachweis der Konvergenz des Verfahrens und der Eindeutigkeit der Lösung 


Aus (3) folgt 


ne | 
Ku) = Mu AARR Ma nen en 


Betrachtet man den ersten Term der rechten Seite von (12) für sich, so sieht man, daß er die 
Lösung der der ursprünglichen Differentialgleichung (1) zugeordneten „verkürzten linearen 
inhomogenen Differentialgleichung 


Be I Be 
x (0) I C; 
im Bildbereich ist. Für (13) ist aber eine eindeutige, rücktransformierbare Lösung gesichert, 
so daß man: 
M,„.(s) + Z(s) 
Up ad Lu Aal a A 
Q | ann 2) (14) 


schreiben kann. 

Aber auch der zweite, rechts stehende Term von (12) läßt sich rücktransformieren. Es ist 
D(s) ein Polynom n-ten Grades in s. Also ist 1/D(s) eine gebrochen rationale Funktion von s, 
und für eine solche ist die Existenz einer Originalfunktion ö(f) gesichert: 


ar A RETTTT HIER 


Ferner hat nach den bisher gemachten Voraussetzungen 2 {f(x,, &,)} eine Abszisse s, > o der 
Konvergenz. Und sowohl ö(t), das die Form 


Ne ; 
' r er et 
PN ze, 2 I a ee 


hat (wenn D(s) I< n Wurzeln o, jeweils von der Vielfachheit m, hat), als auch /(z,, &,), welches 
die Voraussetzung (6) erfüllt, gehören der von G. Dorrsc# [5] definierten Klasse der /,-Funk- 
tionen an. Damit ist aber auf 

h 1 

U u} Die) I u) } 20} 


der Faltungssatz anwendbar: 
s a - x \ . 
Ua u) KÜ = SIEG) z(r)) d (E— 7) dr) a ie A 
0 


und die Rücktransformation gibt: 
[ 
ie KH HD) 


Damit führt die Rücktransformation der gesamten Gleichung (12) auf die schon von L. A. Pıpzs 
verwendete Integralgleichung 


+1) = zit) + [ot Kl); le ER 


die hier als Iterationsvorschrift geschrieben ist. Es ist eine nichtlineare VOLTERRAsche Integral- 
gleichung. Sie ist völlig gleichwertig mit der Iterationsvorschrift (4) des eigentlichen, von P. J. 
NowacktI vorgeschlagenen Verfahrens. Konvergenz von (4), Existenz und Eindeutigkeit einer 
Lösung, die Element des Raumes X ist und die gegebenen Anfangswerte erfüllt, müssen sich wegen 
der eindeutigen Umkehrbarkeit der Lartacr-Transformation also auch an (19) zeigen lassen. 

‚ Auch die Iterationsvorschrift (19) führt für stetiges und stetig differenzierbares x;(!) 
weil ö(l) und x (f) stetig und glatt sind und weil die Integration glättend wirkt, für x(d; +1 nicht 
aus dem Raum der stetigen, glatten und damit dem Raum R der stetigen und stetig difroreniniert 
baren Funktionen heraus. Die Vorschrift (19) ist daher ebenfalls in R unbeschränkt durchführbar, 
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Nun sei zu einer Funktion he #t eine Funktion ge R definiert durch: 


t 
Sy + FEATURE N... 2... (20), 
so daß mit dem Operator 8 gilt: 


gehe NE NIE (21). 


Es ist zu prüfen, ob dieser Operator lipschitzbeschränkt ist. 


Der Operator % führt, wie schon gezeigt wurde, nicht aus R heraus. Es ist R ein metrischer, 
vollständiger Raum, für den man eine Norm 


„Nalul + Nalil 
ert 


lull = Ma ‚ für we, 


N, N, Lirschtrzkonstanten von f, aa la 
) eine passend gewählte positive Konstante 
definieren kann, 
Man erhält aus (20): 
t 
%®— 9 nf LT Te ae (23), 
was mit 
Hal, = N Ne vn ee, (24) 
gibt: 
c 
re A Fe 
oder ausführlicher geschrieben: 
t 
H= A dt — 7) he), Ae))— Me), ahn)d ...:.... (26). 
Damit findet man mit =t-—r: 
U 
dr döfE 
H-- | OO en). 


0 


Weil D(s) mindestens vom 2. Grad in s ist, d. h. also der Differentialoperator in (1) min- 
destens von der 2. Ordnung ist, wird aber ö(0) = 0. Der Beweis hierfür ergibt sich folgender- 
1 


maßen: Aus (16) ersieht man, daß d(0) = 3) ai ist. Die af treten als Zähler bei der Partialbruch- 


r=1 
zerlegung von 1/D(s) auf. Zu der Bestimmung dieser Zähler erhält man eine Reihe von linearen 
Gleichungen, von denen gerade eine lautet: 


1 
2a=d, 
ve 
so daß für n > 2 stets ö(0) = 0 ist. 
Man erhält somit einfacher: 
t 
B dH dö(E) 


0 


und ersieht aus (28) sofort, daß H stetig ist. Dann ist aber He R und man kann auf H ebenfalls 
die Norm (22) anwenden, was später erfolgen soll. 

Es ist 6 (E— r) wegen 0<t,r < T und der bereits nach (16) festgestellten Stetigkeit und 
stetigen Differenzierbarkeit von ö(f) in diesem Intervall samt seiner Ableitung beschränkt, so daß 


lölt, z)| <N;, 
dö te See este a 


gilt. 
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Aus (25) ergibt sich dann wegen (5) und (29): 
! ®. 
IH] sN;,/[IN |, — h,| +N, Ihe — hul] dt ae N Ce: ae (30), 
ö 


und unter Verwendung von (22): 
t 


N 
IH| < x. [e dt ||, — h,|| = her — 1) —ullı #502 0% (31). 
ö + 
Ganz entsprechend ist gemäß (28), unter Verwendung von (5) und (29): 


5 t . . 
IH] s N, f [N, |, — h,| nn N, [hg — h,| dt ns  . er Dee (32), 
0 


woraus mit (22) folgt: 
1 


. N 
A| < N. [ala Ze Dh . HER FRE 
0 
Wendet man (22) für H an, so erhält man mit (31) und (33): 


NN; + NN, 
. 


Aus (34) ersieht man, daß der Operator $ tatsächlich lipschitzbeschränkt ist mit der LiPscHTrz- 
konstanten: 


IIa]| < TR © 


na lo > rer (35) 
und daß 
- KRr,e1l für AND; + NN Er 


ist. Da / frei wählbar ist, kann man die letztgenannte Bedingung für A immer erfüllen. 
Man kann daher unter Verwendung eines Satzes von L. CoLLaTz [6] schließen: 


Unter den Voraussetzungen a)—d) stellt die Iterationsvorschrift (4) oder die ihr gleich- 
wertige Iterationsvorschrift (19) eine eindeutig erklärte Operation in einem vollständigen Raum R 
mit einer LırscHırzkonstanten K, <1 dar. Dann hat aber die im Raum Rt unbeschränkt 
ausführbare Iteration eine eindeutige Lösung, die auch aus dem Raum X ist. Für eine Fehler- 
abschätzung steht die Beziehung 


E—u+ll 2 — — kr +ı — x;|| ee ee ne (37) 
L 


zur Verfügung. 


6. Fehlerabsehätzung. Ein Beispiel 
Die Fehlerabschätzung ist grundsätzlich nach (37) möglich, doch läßt sie sich im Einzelfall 
oft verbessern. Das soll jetzt an einem Beispiel gezeigt werden. 
Voraussetzung: /(x,, &,) sei von der besonderen Form: 


a) = PIE HP) or nr nn 0 en 5 
die Differentialgleichung sei mindestens von 2. Ordnung; es existiere 
[.oi2,) dm, = De Er a ae RE 


D(x,) und y(x,) seien LArLAcr-transformierbar als Funktionen von { und lipschitzbeschränkt 
als Funktionen von x: 


Da) — Da) ZN; ra — al | 
Bad) — ya) <SN a — al J 


Die letzte Annahme ist möglich, weil sich alles wieder in dem abgeschlossenen Bereich 
® : 0=s IST;,A<Sx,&, SB abspielt und man beliebige, aber passende Festlegungen über 
die Fortsetzung von ® und y außerhalb von B machen kann. 


ii! 
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Dann folgt durch partielle Integration aus (19): 


t f 


241) = a) + BL) Sa) + | DO) Hay) dr + | KUH) ylalı))dr (1), 
| 


ds 
ö 
mit } €E=t1—1r. 
Aus (41) erhält man: 
t t 
% +1 % = he [da — Dar ı)] dr + IE dt — 7) [y@&) — van ı)l dr. . (43). 


0 
Daraus folgt mit (29) und (40): 
t 
mal S(NN A NN) u a-ıld N A Mo. 5 (44). 


Einem Vorschlag von J. SCHRÖDER [7] folgend, sei der Operator 


DE ee dr RE Re (45) 
eingeführt. Mit der Abkürzung 
NS NINEN, Nee ar ne ee ee (46) 
schreibt sich dann (44): 
is Del N Pole al Aa (47). . 


Mit (47) hat man eine bessere Möglichkeit zur Fehlerabschätzung, denn, wie man leicht 
einsieht, ist mittels wiederholter Integration: 


- k—2 
I em 
|& +1 — &,| S Nr-1 sn Bo 2 dr, (Kr 2, 2. 2.100). la (48) 
0 
und damit 
t 
«(d) — | SN / ee ln) Dlder (49). 


Es ist ein Ausdruck, der den Fehler zwischen der wirklichen Lösung x(f) und der zweiten 
Näherungslösung z,(f) abschätzt, wobei x(!)e # die Ausgangsstufe der Iteration war. Da der 
tatsächliche Fehlerbetrag mit wachsendem £ immer größer und die Abschätzung (49) immer 
gröber wird, muß man für große Intervalle 0O<t< T mit höheren Näherungen x,, (k > 2), 
und mit der Abschätzung 


cd) — x.) = N [ BIC FE LT) et ld) den ee (50) 


arbeiten. 
In dem von P. J. Nowackt [2] angeführten Beispiel war die Differentialgleichung 


t{+2:+7r= —- 10 2r — 102 
gegeben. 
Also sind 
o= 10%, 
B=-—,., 
yv=—-1I. 


Für die Anfangsbedingungen x(0) = 0, &(0) = 1 wurden nach der Vorschrift (4) die Nähe- 
rungslösungen 
sd =te" 
2 = 875 — Hd et — (8,75 +514+ 23,52) e®! 
gefunden. 
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Es ist 


Dis)=#2+2s+l1l, 
also 1 ih 1 
De) W@+ IP 
und damit dd) ee temt, 
Aue) a ee 
5 Le | 


Trägt man die Näherungslösungen graphisch auf (siehe Bild), so sieht man, daß für =) im 
Intervall 0O<t< 0,6 etwa die Schranken A* = 0, B* =0,33 gelten werden. Alles übrige sei 
ebenfalls auf dieses {-Intervall bezogen: 

N, = Max |ö| = 0,33, 

N, = Max |1 tt Netm...|=10 0ST206, 


N, = Max Te 
0x | 


N, = Max = 0-0 = 333 
0x 


Also ist 
N=11-140#033-33=11+109=219: 
Aus dem Bild ersieht man noch, daß 
I, — | < 0,051 


ist. Damit folgt nach (49) für den Fehler an der Stelle 
t= 0,6, der für 0 <t< 0,6 auch der maximale ist: 


0,6 
le — 2,|< 2,19 [ 106-9 0,05 rdr, 
0 


also 

e —%| < 0,03. 
Da x,(0,6) » 0,3 ist, wird x(0,6) < 0,33, also B* » 0,33 sein, so 
daß die Schätzung O<Sr< 0,33 für den x-Bereich richtig war 
und die Rechnung daher nicht wiederholt zu werden braucht. 


7. Zusammenfassung 


Es ist gezeigt worden, daß unter bestimmten Voraussetzungen und festgelegten Anfangs- 
bedingungen die Lösung von nichtlinearen Problemen mittels LarLAackz-Transformation auf 
Grund eines Iterationsverfahrens für ein abgeschlossenes Intervall beliebig genau erhalten werden 
kann. Die Fehlerabschätzung gelingt durch Zurückführung auf eine gleichwertige Integral- 
gleichung. 

Wenn auch nicht Lösung und Fehler für das nach rechts offene {-Intervall angegeben wird, 
was den Praktiker am meisten interessieren dürfte, so ist es doch immerhin möglich, bei Ver- 
wendung höherer Iterationsstufen, für die rechte Grenze des {-Intervalls genügend große {-Werte 
zuzulassen. Man wird daher in der Lage sein, zumindest die Überregelung auch im nichtlinearen 
Fall sehr genau abzuschätzen, wenn man sich über die Stabilität des Vorganges Kenntnis ver- 
schafft hat. Diese Möglichkeit dürfte für die Praxis nicht ohne Bedeutung sein. 
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‚KLEINE MITTEILUNGEN 


Die Riffelung eines von Wasserwellen über- 
spülten Sandbodens _ s 


Im Folgenden soll eine einfache Erklärung dafür 
versucht werden, daß der Sandboden des Meeres im 
küstennahen Flachwasser oft ziemlich regelmäßige 
Riffelungen aufweist. Es soll nämlich gezeigt werden, 
daß unter einfachen und plausiblen Annahmen eine 
leichte Welligkeit des Bodens instabil ist und mit der 
Zeit anwächst. Hierbei wird vorausgesetzt, daß das 
Wasserin Bodennähe in ständig strömender Bewegung 
ist und daß die Geschwindigkeit periodisch mit der 
Zeit ihr Vorzeichen wechselt. Die zeitliche Periode 
sei aber so groß, daß man die Vorgänge in Boden- 
nähe als quasistationär betrachten kann. Eine 
solche periodisch wechselnde Geschwindigkeit in Bo- 
dennähe wird durch die langen Dünungswellen er- 
zeugt, deren Wellenlänge von derselben Größen- 
ordnung oder größer als die Wassertiefe ist; eine 
weitere Ursache für periodisch wechselnde Boden- 
geschwindigkeiten sind die Gezeitenströme. Gleich- 
zeitig ist die Wellenlänge bei Dünungs- und Gezeiten- 
wellen sehr viel größer als diejenige der Boden- 
riffelung, eine Voraussetzung, die im Folgenden im- 
plizit gebraucht wird, weil das Geschwindigkeits- 
profil in Bodennähe als unabhängig vom Ort an- 
genommen wird. 

Der Transport von Sandkörnern muß durch Form- 
und Reibungswiderstand der Körner im bodennahen 
Geschwindigkeitsfeld bewerkstelligt werden. Wir 
machen die plausible Annahme, daß die Zahl der in 
Strömungsrichtung transportierten Sandkörner pro 
Zeiteinheit und Längeneinheit quer zur Strömungs- 
richtung oberhalb einer gewissen kritischen Wand- 
schubspannung linear von der Wandschubspannung 
abhängt; unterhalb der kritischen Wandschub- 
spannung seien die Körner in Ruhe. Außerdem 
nehmen wir an, daß diese kritische Wandschub- 
spannung während des periodisch wechselnden 
Strömungsvorganges überschritten wird, damit über- 
haupt ein Transport von Sandkörnern auf dem Boden 


stattfindet. Änderungen der Wandschubspannung in 


Strömungsrichtung entsprechen dann Anhäufungen 
oder Abtragungen von Sandkörnern, und zwar hat 
man bei 2 > 0 Abtragung und = < 0 Anhäufung 
(t = Wandschubspannung, & = Koordinate in Strö- 
mungsrichtung). 

Wir nehmen nun weiter an, der Boden sei leicht 
gewellt, die Bodenform sei durch 


Fe en ee VE N) 


beschrieben (vgl. Bild). 


iz 


Die durch das darüber hinwegfließende Wasser 
verursachte Wandschubspannung besteht dann aus 
zwei Anteilen [1,2]: einem Anteil 7,, der auch ohne 
Wandwelligkeit vorhanden wäre und der von x nicht 
abhängt; außerdem aus einem Anteil ,, der durch 
die Welligkeit verursacht wird. Hierfür kann man 
schreiben [1, 2] 

u =02008 (ke Ha). in.» .. (2) 


a>0$ ist eine von Zähigkeit, Grenzschichtprofil, 
Wellenlänge 2r/k und Außengeschwindigkeit ab- 
hängige Konstante, x eine von denselben Größen 


abhängige Phasenverschiebung des Anteiles r, der 
Wandschubspannung gegenüber der Wanddefor- 
mation. | 2 
Da vorausgesetzt wurde, daß die Außenströmung 
periodisch ihre Richtung wechselt, wird man nach 


einer halben Periode folgende Wandschubspannung - 


haben: einmal den von x und 2, unabhängigen Anteil 
— 7, und den von x und 2, abhängigen Anteil 77: 


y=—anca(hr—d) ....6). 


Im zeitlichen Mittel über viele Perioden hebt sich 
der Anteil 7, aus der Wandschubspannung heraus 
und es bleibt nur ein von r, und r/ herrührender An- 
teil, für den man ohne ins Detail zu gehen ansetzen 
kann (— bedeutet zeitliches Mittel): 
er rru+m. 
elle): 

Infolge von Anhäufung oder Abtragung von Sand- 
körnern wird z, sich langsam mit der Zeit ändern. 
Nach dem, was oben über Anhäufung von Sand- 
körnern gesagt wurde, hat man für die zeitliche Än- 
derung von 2 anzusetzen: 


tm —az,sinkxsina . 


dz dt 
en TEE AO), 
aba, 1 
cos ka 7° am coskasin a 
oder ’ 
20 i 
Tann nn (6). 


Offenbar wächst die Amplitude-z, der Bodenwellig- 
keit mit der Zeit an, wenn 0 <« << 180° ist. 

Über die Größe von « liegen verschiedene Unter- 
suchungen vor [1,2]. Im einfachsten Fall einer 
laminaren Strömung, deren Geschwindigkeit linear 
mit dem Abstand vom Boden anwächst und unter 
der Voraussetzung k2, <1, liegt « je nach Größe 
der dimensionslosen REeynoLps-Zahl R = (dU/dz)/v k? 
(dU/dz< = Geschwindigkeitsgradient, v = kinemati- 
sche Zähigkeit des Wassers) im Bereich 0 <« < 30° [2]. 
Für diesen Fall kennt man auch die Abhängigkeit 
der in (2) eingeführten Größe «a von « und kann daher 
denjenigen Wert von R ermitteln, für den das Pro- 
dukt a-sin« und damit die Anfachung der Boden- 
wellen am größten ist. Man erhält für das zur maxi- 
malen Anfachung gehörige R einen Wert von der 
Größenordnung 1. Um hiermit Riffelungen von 
einigen cm Länge zu erklären, wie sie im küstennahen 
Flachwassergebiet oft beobachtet werden, muß man 
bei der Annahme, daß dort dU/dz von der Größenord- 
nung 1 sec! ist, einen Wert von v * 0,1 bis 1 cm?/sec 
einsetzen. Ein so hoher Wert der kinematischen 
Zähigkeit könnte in Wasser allerdings nur durch 
turbulente Vorgänge erzeugt werden. 

Diese Überlegung hat natürlich keinen quanti- 
tativen Charakter, z. B. deshalb nicht, weil das Ge- 
schwindigkeitsprofil sicher nicht als linear voraus- 
gesetzt werden darf. Man kann aber auch für ein 
allgemeineres laminares Geschwindigkeitsprofil zeigen, 
daß «x im Bereich 0 < «a < 180° liegt (vgl. z. B. [1]). 
Andrerseits ist zu bedenken, daß in Wirklichkeit 
turbulente Vorgänge die Strömung in der Bodengrenz- 
schicht bestimmen werden. Aber auch bei turbulenten 
Grenzschichten ist es jedenfalls plausibel, daß x in dem 
Bereich liegt, der Anfachung ergibt. 

Schließlich sei erwähnt, daß es verschiedene andere 
Deutungen der Riffelung eines Sandbodens gibt [3, 
4, 5]. Eine ausführliche Übersicht findet man in [3] 
und [4]. U. U. kann die oben erwähnte Instabilität 
die Ursache für die Bildung von Sekundärriffeln im 
Sinne von [4] sein. 
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Stresses in a Twisted Cylindrical Plug of Iso- 
tropic Material Inserted into a Rigid Medium 


Introduction 


In this paper we have considered an elastie plug 
in the form of a short eireular eylinder with one end 
and the lateral surface held fixed while the other 
terminal section is acted on by shearing forces dis- 
tributed over a small circular area. Two different 
types of distribution of shearing forces over the eir- 
cular area namely (i) linear distribution and (ii) con- 
stant distribution have been discussed. Numerical 
caleulations of shearing stresses at the fixed end have 
been computed for different values of r, the radius of 
the plug. 


Formulation of the problem 


Let r, 06, z be the cylindrical coordinates with the’ 


axis of the cylinder along the z-axis. The origin being 
at the fixed section and assuming = u=0(0 and 
ug = v to be independent of 9, we find that the strain 
components are given by 
= 
and the stresses are 
"=-b=-Rr-r=0, 


Tg EL 2 Du 2). 
a=u[, r). az Role 


The boundary conditions for the eylindrical plug 
of length h and of radius « twisted by means of a 
shear distributed over a part of the section z= h 
with the end 2=0 and the lateral surface r= a 
held fixed, are given by 


Hemden m —0, 


ov ® ov 
ee ER TE 


&6 = — 
i ar 02 


Zell, er 

0s2z sh, v=0, ‘ zu9)« 
ht byerreares Kr) 

From (2) it is found that two of the equations of 
equilibrium are identically satisfied and the third is 
given by 


VERds, 


02V 1 % v 0% 
or: r or ry2 


R m (baue 
Solution 


It we goin for product solutions of the form Z(z) K(r) 
which satisfy the boundary condition 


YW,o=0, 
one type of such solutions can be written as 
v= ERılr)sinh (k2) . 2... .6), 


where k is an arbitrary constant and R satifies the 
differential equation 

d?Rr 1 dRr B 1 
ee (r— =) u (6). 


The solution of (6) is 
Bil) = Jen) » 0... 00% 
where J, is the BzsseL function of the first kind and 
first order. ö 
"Therefore from (5) and (6), we get 


v=AyJılkr)sinh (k2) . . . - (8). 


Now from the boundary conditions (3), since atr = a, 
v= (0, we get 
Aka) 0 RT ar 


Let kysky kp... > kn be the roots of the equation (9). 


Then 


an & Andılkar) Ai (ke)... Zr 
and ” 
(Den = (r BE), kn An Jılknr) cosh (ku) 


=h 


= Fr), for, 0 sr Sb... (1% 


so that 
se Bu) 
” 7 pe u kn cosh (kn h) [I (km b)]% 
E 2 
x [ Fan nn de Me (12), 
0 
b 
=B, hi r Fir) Jılkar) dr, (12.1) 
r | 
where 
2 
= (12.2). 


B? 1 kn cosh (kn A) [I 4 (km DE 


We are next to consider several possible values of 
Fr) depending on the different modes of distribution 
of shearing tractions. 


(i) Linear distribution of shearing forces over the 
terminal section not held fixed. 


In this case 
Ar)= 8% Vzersh, 2 un 
=. Qi; b zuh, 


S being a constant. Substitution of (13) in (12.1) 
yields 


(13). 


b 
An TuS 2 Jlkar)dr = 78 


0 


where J, is the Besser funetion of the first kind and 
second order. 


From (11), (12.2), (14) we get 


on oo 
92 = u 2 kn AnJılkn r) cosh (ky 2) 
Nn= 


© J(knb) Jı(kn r) cosh (kn 2) 


=2$8-7 . sh (‚ R 
n=ı An cosh (k„ h) [J% (kn b)]? (15). 
From (2) and (10), we have 
be er hi 
For > A„sinh (k, 2) ji Iilln ) = re (16). 
“ n=1 r 


Substitution of A, from (14) and (12.2) in (16) yields 
SEELEN 
k3, [J1(kn b)] 


n=1 
’ J(knr)) sinh (ky z) 
| nl) r  jeosh(knh) 


BEN ) \ Jz(kn b) Jy(kn r) sinh (Ay z) 


(17), 
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and , 
UN En S Jz(lkn b) Jz(kn a) sinh (ky 2) 
a oe 


n=1 


(ii) Constant distribution over a circular area. 
Here 


(18). 


Ares 0er», 
0b ng 


llin 
z=.h, 


where Sisa constant. Substitution of thisin Eg. (12.1) 
gives 


EA) 


b 
An= TnS[r Ilka r) dr. 
0 


The integration can be carried out by expanding the 
B#sseL function into a series of ascending powers 
of k„r and integrating the result term by term. 

Thus 


b 
+ AR Kuren +1 
er) 2 
0 
or 
[eo] > 
}o (— 1)r k,\2n+1 b2n+3 = 
er no nach 2) 2n +3 (20), 


where 7'„is given by Eq. (12.2). 
The stress components are then given by 


an oo 
9 = u DD knAnJılknr) cosh (kz) . . (21), 
n=1 
and 


2 © 
(r0),-a = — 4 3 km An sinh (kn 2) Ja (km a) (22), 
n=1 
where A„is given by Eq. (20). 


Numerical caleulations 


Numerical values for stresses 62 in the first case 
have been obtained at the fixed end 2=(, for 
different values of the distance r. The numerical data 


considered are «= 1,b = and h=2. Variation 


of the distance r thus lies between 0 and 1. Results 
obtained as shown graphically in the accompanying 
diagram indicate the stress (92)z-0o inereasing with 
distance from the fixed axis and reaching a maximum 
value at r=a=1|, i.e., at the periphery of the 
isotropic plug. 


1 ai 
0,50 075 


Lr—- 


) an 
0 0,25 


Showing stresses at the fixed end forlinear distribution 
of shearing tractions at the other end 
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Note on the Bending ofa Thin Uniformly Com- 
pressed Circular Plate under an Eccentric Con- 
centrated Load 


In a previous paper (BAsuLı, 1960) the problem of 
bending of thin clamped circular plates submitted to 
the simultaneous action of uniform compression in 
the middle plane of the plate and concentrated load 
at the centre or any other type of load which is a 
function of r has been solved. Here the corresponding 
problem with an eccentrie concentrated load is being 
discussed, the boundary of the plate being clamped 
in this case also. The solution is based on the approxi- 
mate plate theory. 

Except at the load, the deflection w of the plate 
satisfies the equation 


WE Re) ee (1), 
where 
02 02 
ar u ES 
vi ox2 a: oy2 
ad E R 


— Compression per unit 
length of the eircular boundary, 
2Eh: 
31 ae) 
E modulus of elasticity, 
2h = thickness of the plate, 
o = Poıssonx’s ratio. 


D= ; flexural rigidity, 


I 


lt N, be the shearing force on a circular ring of radiusr 
then N, is given by 


ö 2 < 
Eu e a (2): 
m 
For clamped edge, we have 
ow 
el a RN) 


on the boundary. 

Let the concentrated load P be placed at a distance 
b from the centre of the plate and the radius of the 
plate be a. To solve the problem let us divide the 
plate by a concentric cylindrical surface of radius 5 
passing through the load. Taking the line joining 
the centre of the plate and the load as the initial line 
and centre of the plate as pole the deflection equa- 
tion (1) can be written as 


02 1650 NaoeH 
(+ ” 


r or r? 002 


I at u)=0 ER: 


or rar "nr 00 


we ze 


As solutions of (4) we assume 


oo 
w=w=#%+ & Rmecosm6, whenr > Deal), 
m=l 
and 
= e) 
wv=ww=R,+ 3 Rmcosm6; when r<b . (6), 
ml 


where Ry Ri, Km and Ri, are functions of r only. 
15 
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Substituting equations (5) and (6) in (4) anp con- 
sidering the contributions of R, and R% alone we shall 
see that they satisfy an equation of the form 


2 

(Be run) m 

Solutions of this equation may be put as 

R= R, = A Jılar) + Bo Yılar) +0, + Dlogr, 
Eichen: 72 DE Se LE (8), 


R=R=AyyJılar) +05, whenr<b...(9) 


Jı(&r), Yo(@ r) being Brsskr’s functions of dar zero. 
If the outer boundary be clamped we have 
ow, 


w=—=0, 


Er when n=a = ...(10) 


As the deflection, slope and bending moment will 
be continuous on the dividing eircle we get the 
relations 
ow, dw, , rw, _ ÖFw, 


ee = a whenr=b(ll). 


BR 


Equations (8) and (9) contain altogether six con- 
stants and relations (10) and (11) are five in number. 
To get the sixth relation we shall have to consider 
the shearing force on the dividing ceircle which is 
continuous at every point of that circle except at the 
concentrated load. Representing the load in the form 
of an infinite series 


e + e Ka) 


the discontinuity in the ua force is given by 


— D>- 
r=b 


0) 
pefosrejm| otutsene 


ee > u 
2 = 


Substituting (8) and (9) in (10), (11) and (12) and 


solving for A,, By, Oy, Do, Ay and 0% 
122 
er T xaJ,(«a) 


j a Ju(& b) [Yı(x a) Jular) — Yolar) Ja a)] 


—J,(za) —xaJ,(xa)log 3), 


+ Joan) + Julab) 


1 


I - 
2x TaaJ,(«a) 


0 
l a Jar) [Yıl@ a) Ja 6) — Yolz b) Jıla a)] 


+Jolar) + Jlab) —Jol@a) —xaJ,(xa) log , K 

) 
Substituting equations (5), (6) in (4) we find that 
equations to determine Rand Ri) will be of the form 


Kl ” N ° a m? 0&@R l öR ae r 
or Or  va)\ Ort h rfr pi A: R) er 
RT We 18), 
This equation is satisfied if 
R= Rm= AmJm(&r) + Bm Ynlar) + Cmr® + Drn; 
whenirge Den Sn, ae (16) 
R= Rn = AyJmla N Omr® whenr<b, . (12) 


where "Jn(a, r), 


Ym(xr) are Besser’s functions of 
order m. 


Substituting a and (17) in (10), (11) and (12) and 
solving for An» m; Din, Am and O7, we obtain 


r 


a T aa Jm+ı(“ a) 


— Yan) Imlaa)l+ (ame n) 


+ nn" 


Bu [ra/nta Ya a) Imler) 


v pP 
Rn ra Talea) |? Imla r) LYmsılza) Inka 6) 


a le ae 21 Imlar) 


Y Ant = (ef ge Imılaa) 


__[a\mxa ]/ r\® 
) 
The deflection w, of the outer portion can be ob- 
tained by combining (5), (13) and (18) and w, that 
of the inner portion can be obtained from (6), (14) 
and (19). 
P 
x 
2n Txal,laa) 


Central deflection = 
l al Y (a a) Ju(@ b) — Yılz b) Jı(z a)] 
+1+J,(&b) — Jıl@a) — za J,(za)log ’) (20). 


When the compression vanishes we have on mak- 
ing & tend to zero, 


0 2 HE) 
+ 5?) log > 


20°, 


Dr b2) log 


b ze + al u a 


Pb 
16x D 


BATETERT er 
at b2 9, 


Pb» (a? — 2)? 73 


I, en ji 


a? b? ab 
4 E [77 
’ 


 1WrD Trade 
E87 


Vi hm pm 
3 m(m — 1) 7) |" mn 1) 


Rn= 


m(m — 1) b? | 
m+1 =) 


I m—|1 
Tom ( Kat ) ’ 


Y P pm yMm 
”" Smm—1)rD \omr 


arm PM+2 


217 m ” +(m—.1) 


+ (m l)r?- 


(m — 1) 


(% Fr 


m 1 
m m + 2 _ m = y"ı 


s 
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Combining these the deflection at the load will be 


P (a — 522 
16r.D a? 


' "This result tallies with that of Tımosuzxko (1940, 
page 268) obtained in the case of a plate without 
compression. The central deflection in this particular 
case reduces to 


IE b 
2 Br 
ll An 2 


The deflection w of the plate under the simultaneous 
action of uniform compression in the middle plane of 
the plate and a concentrated load at the centre of the 
plate can also be deduced from (5), (13) and (18), 
ra limits as b tends to zero, in the form (BAsULı, 


2 
Nas). 


BR 
W= DET aa Ilaa) l al Yı(@a) Jar) — Yılar)Jı(@a)] 
+Jo(@ r) +1 — Jo(la a) — a a J,(a a) log 2 (24). 


The corresponding central deflection is 


P 
2rTaayıaa) Aa) 


IR EInle a la a)logie a) sh. 


As expected this result reduces to 


Pa: 


Br 9 
I6xD when &—>0 ..... (26). 


In conclusion author wishes to thank Prof. B. Sen, 
D.Se., F.N.I. of Jadavpur University for his gui- 
dance in preparation of this note. 
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Flow of Visco-Elastic.Maxwell Fluid through 
Tubes under Exponential Pressure Gradient 


Introduction 


Flows of NewTonian viscous fluid through circular 
and elliptic tubes have been discussed by several 
authors [1]. Here we have considered the flow of a 
visco-elastic fluid of the MAxwELL type (i. e. a spring 
and a dash-pot arranged in series) governed by the 
equations 


W;=—Ppsit ti 
BI, 
Hr a) = un, (1), 
1 
;=5liit vi) 


where 7;,;is the deviatoric stress tensor, p is the pres- 
sure, x, u are material constants (x being of the 
dimension of time is called the relaxation time und u 
is the viscous parameter) and v; * = 1,2, 3) are the 
velocity components. If x is large the fluid behaves 
more like an elastic solid and if x is small the behaviour 
is like a viscous fluid. Two types of laminar flow will 
be discussed, 1) flow through a;circular eylinder, 
2) flow through an elliptic cylinder. 
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I. Circular Cylinder 


Choosing cylindrical polar coordinates (r, 0,2) with 
2-axis along the axis of the cylinder, the components 
of velocity are 


‚a=0, v=0, w=ufrt) 


(2): 
as the motion is assumed to be laminar, and the 
equations of motion give with the help of (1), 


8\ dw 1 8\ m 
DH a) ne ie 
| le ee 


u 1m 
Fe „le () 
1 9» 


The other equation vanishing identically from sym- 
metry. By virtue of the second equation of (3) we can 
assume 


In particular we take 
lt) =ae-mit, 
where «, m are real constants. In that case we assume 
tentatively 
w = ner Wi ra ne a (8) 
with the boundary condition 


f(a) = 0 

Equation (3) then simplifies into 
d? ld 

Prag 


m? (1 — x m? 
dr? r dr v 


") fr) + am] =0 . (). 
Solution of (5) is, using (4a), 
fl ala] ee ii 


where J, is the BEssEL function of the first kind of 
zero order. 

We now consider two cases: 

Case a): E |na| <1, then using 


= 
Ja) al — TE 


when x is small, we get from (6) 


a(l — x m?) 


Me EN: 


This result differs from that for an ordinary viscous 
flow in the factor «(1 — x m2)/4»v in place of «/4v 
showing that the presence of the constant x flattens 
the velocity profile. This is expected as the elasticity 
of the fluid damps the flow. 

Case b): I x is large such that «m? >1 and 
n'a>1 (writing —n’? for n? in (6)), 


N = | | FRE NSN 


I(n’ a) 


where I, isthe modified BzsseL function of the first 
kind of zero order. 


But as 
I,(x) — exp (&) Ä +o Bl Are), 
when x is large, so (8) gives 
fir) = 5 (alr!* [exp {m ta} 1]. . O), 


showing the flow has a boundary layer character. 
15* 
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Elliptie Cylinder 

We take z-axis along the axis of the cylinder and 
x, y-axes along the major and minor axes of the 
ellipse. In that case 


vw=0, v=0, w=ui,yP). 
Then assuming as before 
_ 1%» =NX e-m*t 9 
0 82 

the equations of motion give 

&\ dw : rw , uw 
( +e) F7 — «(l — x m?) em t a +) 

(10). 


Assuming w(z, y, t) = f(x, y) e-m*t with w= 0 on the 
boundary (x?/a?) + (y?/b?) = 1 equation (10) gives 


Fi 


en m? (L— x m?) 
da? " ay? v 


—/remj]=0 (U) 


Changing to elliptic coordinates &,n by the substi- 
tution 


z+iy=ccosh(£+in), c= Ja —b2 


and separating the variables, 


d29 m? (1 — x m?) en 
7a Ar aa cosh 22|0=0 (12), 
d’y m? (1 — x m?) £: 

| en 3 2cos2n]|y=0 (13), 


where we have taken 


= 1+ am? = OS) yln). 
Case a): Assuming, 


al 
a in. 

2v 
and since wis symmetrical about the major and minor 
axes of the ellipse and is single-valued and periodie 
in n with period x, so y is a multiple of the MATHIEU 
functions cea„(n, 9), so that we have 

00 
Ku 2 Con Cezn (8, q) cezn (N, N; 


Nn= 


oo 
ln(,)=L% arm cosh2rE&, 
r=0 


(14). 


oo 
cean (N, 9) = 2Ayr cos2rn 


Using the boundary condition 
zw=em, E=5 

and employing the usual normalisation relations 
; Ogn = a A Im Oegn (Eu 9). 

We know that AP 1, as g—0. So in view of our 
assumption if g <1 we expand [2] Ce,„ and ceyn in 
powers of q and retain only first power. Thus, since 

1 : ’ 
AP=1, AP=7g+0@), Am = o(gt), 


therefore, 


1 
X = — h 5 5 1 (vosh 26 — 0002 7 — cosh 28 


2 
m“ y 


cosh 28 c082 7 
cosh 2&, | 
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so that 
1 (oh 250 — 00027 — oh 24 
v 


4 oh 


showing that here also the presence of x damps the 
motion. a. 

Case b): If x is large such that «m? >1 then 
replacing 29 by — 2g', q’ > 0, we have, proceeding 
exactly as before. ag 


TB ent 


= &0om Oean (&; — ') een (m, —Q)> 
en —gd)=(— 1 (- ır AZ cos2rn, (16). 


Con —g)=(—1" z — pr AEM gosh 2r& 


Using the boundary condition and the normalisation 
relations, and remembering that A = 1, Aa® is 
small when n > 1 we have the asymptotie formula [2] 


2 1/2 
Ce (&, —gq’) » (Faake E ;) %, cosh 


ı 


? Yq’ cosh &— tanh-1 [ion (7 _ 2) > 


nd) Jap aa mn, 


when >1. 
Therefore 


& im 1 ‚\va 
= ze 57 :0/cosh 3 :) 
exp [— 24’ (cosh &, — cosh &)] ce, (7, —Q') » 
giving | 
Mm 


1= zal(eosh Zöofeosh 58) "xp (— 2 Ya (och 


— cosh :) ce, (7, —q) — ı AR NT, 


‚showing that in this case alsothe motion hasa bound- 


ary layer character. 


In eonelusion I express my deep gratitude to Prof. 
B. Sen, D. Sec., F. N. I. for his kind help. 
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Zur Abschätzung der Fortpflanzung der Daten- 
fehler bei linearen Gleichungssystemen nach 
der Formel von Wittmeyer-Collatz 


Bei linearen Gleichungssystemen sind in der Praxis 
die Daten meist mit Fehlern behaftet, die sich in die 
Lösung des Gleichungssystems fortpflanzen. Zur Ab- 
schätzung dieser Fehlerfortpflanzung wird eine Formel 
bewiesen. Es wird gezeigt, daß zwischen dieser Formel 
und einer von WITTMEYER und COLLATZ angegebenen 
Abschätzung ein Zusammenhang besteht. Da in die 
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 Abschätzungsformel die Norm der Inversen der 
'Koeffizientenmatrix eingeht, erweist es sich als not- 
wendig, noch eine Abschätzung für diese Norm anzu- 
geben, wozu ein Iterationsverfahren von ScHULZ 
herangezogen wird. Zugleich wird damit die Frage 
beantwortet, wie der Einfluß der bei einem Elimi- 
nationsprozeß zur Lösung eines linearen Gleichungs- 
systems auftretenden Rundungsfehler abgeschätzt 
werden kann. } 


1. Problemstellung 


Von Seiten der Praxis wird häufig.die Aufgabe ge- 
stellt, ein lineares Gleichungssystem der Form 


Ele 


zu lösen, bei dem die Matrix © und der Spalten- 
vektor c zahlenmäßig gegeben sind, während für die 
Elemente öc;; und öc; der Fehlermatrix öC bzw. des 
Fehlervektors öc nur Schranken Ac;; bzw. Ac; für ihre 
Beträge bekannt sind 


da Aa lock 2 Acıl.. 2.2): 
Numerisch können wir nur das Gleichungssystem 
Ge 0 ee (3) 


lösen, dessen Matrix C' als nicht singulär vorausgesetzt 
wird. Für das System (3) sei eine Näherungslösung & 
mit dem zugehörigen Defekt r bekannt. Dann genü- 
gen Zund r der Gleichung 


08 —c=—r 


RIER (4). 


Gesucht wird eine Abschätzung für die Differenz 
zwischen der exakten Lösung & des Gleichungs- 
systems (l) und der Näherungslösung 2 des Sy- 
stems (3). 


2. Abschätzung der Fehlerfortpflanzung 


Zur Abschätzung der Fehlerfortpflanzung gehen 
wir von den Gleichungen (1), (4) aus, die uns sofort zu 
den Beziehungen 

((—6ö0)E —-CZi=öhHr, 

(C—- SO) E —-D)=sICEHÖÜiHr... (Ö) 
führen. Für einen Spaltenvektor x = (x!) erklären 
wir nun eine Norm ||<|| und mit ||C|| bezeichnen wir 
die zugehörige Norm einer quadratischen Matrix 


© = (cj). Im Zahlenbeispiel werden wir unten die 
Zeilensummennorm verwenden; dann gilt: 


Ill = max «| ,” IC = max Z jez] . 
® 7 


Aber auch die Verwendung einer anderen Norm- 
definition kann gelegentlich von Nutzen sein. 


Gehen wir nun in (5) zur Norm über, so können 
wir nach bekannten Normgesetzen schreiben: 


I —zl| S |I(C — 8Oy" 1] AlSCI - |Iell + Ildell + In 


Unter der Voraussetzung 
el sell <1 - (7), 
die beinhaltet, daß die Matrix O’nicht fast singulär ist 
und die Koeffizienten des linearen Gleichungssystems 


nicht mit zu großen Fehlern behaftet sind, können wir 
die Norm der in (6) auftretenden Inversen abschätzen: 


II(C — 80) = |IIC(E — EC" SO)T|; 
= INNE 50)". 
Kal 
= = 
1—je=|} - |18C| 
Damit ergibt sich aus (6) sofort die gesuchte Ab- 
schätzung 


en, 


| ISCH - J1zll + Iöell + IIrll 


I: — 21 s 10 


1 — 110]: 160] 


Da zwischen (8) und einer Abschätzung von WırT- 
MEYER und CoLLATZ ein Zusammenhang besteht, 
wollen wir (8) die Formel von WITTMEYER-COLLATZ 
nennen. : 


In der Abschätzung (8) treten die Normen der 
Fehlerglieder öC und öc auf, die wir aber als bekannt 
ansehen können. Denn wegen (2) gilt ja z. B. bei Ver- 
wendung der Zeilensummennorm 


||öC|| = max 3 |öcır| < max 3 Acır,, 
i & do ik 

|Ide]| = max Ac;. 
% 


Weiterhin geht in die bewiesene Formel die Norm 
der Matrix C-lein. Das stellt uns vor eine neue Auf- 
gabe, nämlich Abschätzungen für die Norm der In- 
versen einer Matrix anzugeben. Mit dieser Aufgabe 
werden wir uns unten näher befassen. 


Die Formel von WITTMEYER-COLLATZ beantwortet 
zugleich die Frage, wie der Einfluß der Rundungs- 
fehler abgeschätzt werden kann. Dazu haben wir uns 
die Daten der Aufgabe als genau gegeben zu denken, 
also die Elemente der Fehlerglieder gleich Null zu 
setzen. Dann ist aber & gleich der exakten Lösung x 
des Gleichungssystems (3) und ||öC]|| = ||öc|| = 0. 
Voraussetzung (7) ist stets erfüllt, und (8) führt zu der 
Abschätzung ? 


ie — zil = ]I0 "| - IIrll 


für die Differenz zwischen der exakten Lösung x und 
der Näherungslösung & des Gleichungssystems (3). 


3. Zusammenhang mit einer Formel von WITT- 
MEYER- COLLATZ 


Für das hier untersuchte Problem hat bereits 1936 
WITTMEYER [1] eine Abschätzung angegeben, die von 
CoLLATz [2] modifiziert und verbessert wurde. Cor- 
LATZ beweist die Formel 


—— . (9) 
ce — |0| 
die unter der Voraussetzung 
ey -PRa >02 (10) 


gültig ist. Dabei wird für einen Vektor x = (x!) der 
Betrag |x| durch die Gleichung 


kI= VER? 


definiert, was gerade der Quadratsummennorm für den 
Vektor x entspricht. In (9) treten weiterhin der obere 


Betrag |C| und der untere Betrag |C| einer Matrix O 
auf, die so definiert werden, daß für einen beliebigen 
Vektor & die Ungleichung 


|. el = |0zls Ol. el 


besteht. Nun gilt aber 
2] = le" Oel = ||C"]| - ||C ||, 


(ul). 

IC =] s|lell- el! 
Verstehen wir (11) im Sinne der Quadratsummen- 
norm, so können wir |C| = ||O||, |E| = 1/||0”!|| wählen. 


Mit dieser Wahl des oberen und unteren Betrags 
einer Matrix geht aber Abschätzung (9) in (8) und 
ebenso Voraussetzung (10) in (7) über, wenn wir (7), 
(8) noch auf die Quadratsummennorm beziehen. 


Damit hat sich ein Zusammenhang zwischen den 
Abschätzungen (8) und (9) ergeben. Die Abschätzung 
nach (8) besitzt den Vorteil, daß neben der Freiheit in 
der Wahl der Norm die Verwendung der Norm von 
0-1! für den unteren Betrag der Matrix C zu einer 
praktisch brauchbaren Abschätzung für die Fehler- 
fortpflanzung führen wird. 


Be 
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4. Abschätzung der Norm der Inversen 


Die Formel (8) stellt uns vor die Aufgabe, eine Ab- 
. schätzung für die Norm der Matrix 0”! anzugeben. 
Für die folgenden Betrachtungen setzen wir voraus, 
daß eine Näherungsmatrix X für die Inverse 0-1 be- 
kannt ist. Diese Voraussetzung ist mit einer erheb- 
lichen numerischen Zusatzarbeit verbunden, die wir 
aber in Kauf nehmen wollen, um eine Abschätzung 
der Fehlerfortpflanzung zu gewährleisten. 

Zur Abschätzung der Norm der Inversen ziehen wir 
ein Iterationsverfahren heran, das von Sc#urz [3] 
angegeben wurde. ScHhurz konstruiert eine Matrizen- 
folge nach der Vorschrift 


a AT EL: 
die gegen die gesuchte Inverse von C konvergiert, wenn 
gel —OZHlı<t: su. (13) 


gilt. Die Voraussetzung (13) beinhaltet, daß die 
Ausgangsmatrix X, eine brauchbare Näherungs- 
matrix für die Inverse O1 ist. Als Ausgangsmatrix 
für den Iterationsprozeß wählen wir die von uns als 


bekannt vorausgesetzte Näherungsmatrix X. Für 
dieses Iterationsverfahren wurden von Dück [4] 
Fehlerabschätzungen angegeben: 


als IK - - (19), 
Hs ll - - 08). 


Die Ungleichungen (14), (15) gestatten uns sofort, 
die Norm von 0! abzuschätzen 


1 St Kl | 0), 


2 
Sl |. - - am. 
Die erhaltenen Abschätzungen sind unter der Vor- 
aussetzung (13) gültig; ihre Anwendung setzt die 
Durchführung eines Iterationsschrittes zur Berech- 
nung von X, voraus, was bezüglich des Arbeitsauf- 
wandes etwa zwei Matrizenmultiplikationen ent- 
spricht. 
Eine besonders einfache Abschätzung erhalten wir 
in der folgenden Weise: Wegen (12) ist 


%=X+X, (E—CK,) 
und damit 
Ill = |I&oll + X] - || — CO X] = (1 +) ||X||; 
was uns mit (17) zu der Abschätzung 


2 Koll 
le Is], NL 


führt. Zur Anwendung von (18) ist es nicht mehr 
notwendig, einen Schritt des Iterationsverfahrens aus- 
zuführen. Es ist nur noch g zu berechnen, was etwa 
eine Matrizenmultiplikation erfordert. 

(16) ist die genaueste der angegebenen Abschätzun- 
gen; am einfachsten, aber zugleich auch am schlechte- 
sten, ist die Abschätzung (18). Welche der Ab- 
schätzungen man in der Praxis verwenden wird, hängt 
wesentlich von der Größe von q ab. 


5. Beispiel 


Die erhaltenen Abschätzungen wollen wir auf ein 
Beispiel anwenden. Wir wählen: 


200 40 20 340 
C=:| 45 150 s) : ‘= [mo - 


107521022100 330 


Es möge |öcjj] < 1, |öcı| < 1 für alle i, j gelten. Dann 
ist ci Bi |ldel| <1. Zur Abschätzung nach (8) 
benötigen wir einen Näherungsvektor 2, zu dem wir 
den Defekt r berechnen 


0,99 T 0,005 
z= 12,02], r = |—0,018|. 
3,01 — 0,011 


Dann ist |lz|| = 3,01, ||r|| = 0,018. Zum Vergleich 
bemerken ya daß für Abe der exakten Lösung 
«|| = 3 gilt. > 

Als Ausgangsmatrix für den ScHurzschen Ite- 
rationsprozeß wählen wir 


ee Me 
200 1000 1000 
En a 2 
% = (500 I000 T00|’ Koll = 7000 
—l —l 11 
1000 1000 1000 
Leicht lassen sich die Werte 
19 
ga= EC X = Ip 
11 11 
Ill = 1000’ I — Zoll = 70000 


bestätigen. 

Der exakte Wert von ||C-!|| und die Ergebnisse der 
Fehlerabschätzungen (16) bis (18) sind in der ersten 
Spalte der folgenden Tabelle zusammengestellt. In 
der zweiten Spalte sind die Ergebnisse der Fehler- 
abschätzung (8) bei Verwendung der jeweiligen Werte 
für ||O-1]| eingetragen. 


Abschätzung von 

| ne | Ie-al | 
exakt " 0,0112 0,117 
nach (16) 0,0113 0,118 
nach (17) | 0,0116 0,121 
nach (18) | 0,0161 0,170 


6. Schlußbemerkungen 


Vor eine ähnliche Aufgabe wie die hier behandelte 
wird man gestellt, wenn man die Inverse einer Matrix 
zu bestimmen hat, deren Elemente mit Fehlern be- 
haftet sind. Wieder erhebt sich die Frage nach einer 
Abschätzung der Fortpflanzung der Fehler in die ge- 


"suchte Inverse. Man kann leicht zu einer zu (8) völlig 


analogen Abschätzung gelangen. Liegt eine Matrix © 
mit überwiegenden Hauptdiagonalelementen vor, so 
kann man andere einfache Fehlerabschätzungen an- 
geben. Alle diese Fragen sowie die bereits mitgeteilten 
Ergebnisse sind ausführlich in [5] dargestellt worden. 
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Constantin Weber und Wilhelm Günther, Tor- 
sionstheorie. 306 S. Berlin 1958. Akademie-Verlag 
Preis geb. 38,— DM. 


Der erstgenannte, durch zahlreiche bedeutende Bei- 
träge zur Festigkeitsforschung und Elastizitätstheorie 
bekannte Verfasser hat zusammen mit einem jüngeren 
Kollegen durch das vorliegende Werk über Torsions- 
theorie eine große Lücke in der Literatur ausgefüllt. 
Obwohl schon viele größere und kleinere Abhand- 
lungen über Verdrehungsfestigkeit vorliegen, fehlte 
es dennoch bisher an einer durchgreifenden und ge- 
schlossenen theoretischen Darstellung, welche auch 
für die vielfältigen in der Technik vorkommenden 


. Probleme eine Antwort zu geben vermag. Ausgehend 


von den elastizitätstheoretischen Grundgleichungen 
des prismatischen Stabes und seinem Spannungs- und 
Verformungsverhalten unter dem Einfluß tordieren- 
der Kräfte werden zunächst jene Lösungen der Diffe- 
rentialgleichung 2. Ordnung angegeben, die sich in 
geschlossener Form darstellen lassen. Weitere Lö- 
sungen lassen sich durch endliche und unendliche 
Potenzreihen der komplexen Querschnittskoordinate 
darstellen (regelmäßige oder unregelmäßige Gerad- 
linien- und Bogenvielecke). Weitere Verfahren (Me- 
thode der GreEENnschen Funktion, Kombinations- 
methode, konforme Abbildung, Spiegelungsmethode) 
werden in den nächsten Abschnitten eingehend sowohl 
allgemein als auch anhand praktischer Beispiele vor- 
geführt. Die zugehörigen speziellen Lösungen werden 
angegeben (Sektor der Kreisfläche, Querschnitt mit 
Loch, Rechteckquerschnitt, symmetrisches Kreuz, 
Kreisquerschnitt mit exzentrischem Kreisloch, Zwei- 
bogenquerschnitt, Kreisquerschnitt mit Einschnitten, 
kleine Löcher und Kerben, Innenecken mit kleinen 
Ausrundungen). Für die Berechnung des sog. Flächen- 
torsionsmomentes, das, mit dem Schubmodul multipli- 
ziert, die Drillsteifigkeit des Stabes kennzeichnet, 
werden mit Hilfe von Variationsbetrachtungen obere 
und untere Grenzwerte angegeben, desgleichen für die 
Torsionsspannung. Das Buch wird ergänzt durch 
zwei im Anhang untergebrachte rein theoretische 
Abschnitte über Elastizitätstheorie und zweidimen- 
sionale Potentialfunktionen, sowie durch ein Lite- 
raturverzeichnis. Das Buch richtet sich nicht nur 
an angewandte und reine Mathematiker, sondern vor- 
nehmlich an theoretisch interessierte Ingenieure. Es 
zeichnet sich neben seiner anschaulichen Darstellungs- 
art vor allem durch eine originelle Art der Entwick- 
lung der Lösungsmethoden aus und wird dankbar 
begrüßt werden. 


München H. NEUBER 


C. Pisot et M. Zamansky, Math&ematiques ge&- 
nerales (Collection universitaire de mathömatiques, 
Tome IV). XXIV + 6488. m. 50 Abb. Paris 1959. 
Dunod Editeurs. Preis geb. 4.500 F. 


Das vorliegende Buch ist eine moderne Einführung 
in die Algebra und Analysis für Mathematikstudenten 
der ersten Semester sowie für mathematisch inter- 
essierte Studenten anderer Fachrichtungen. 

Der erste Teil bringt einige Begriffe und Symbole 
der mathematischen Logik und der Mengenlehre, den 
Funktionsbegriffund den Begriff der binären Relation. 
Die logischen Symbole sind im Folgenden jedoch nicht 
mehr als nötig zur Darstellung herangezogen, so daß 
auch der Anfänger alle Teile des Buches ohne Schwie- 
rigkeiten liest. 

Der zweite Teil ist der Algebra gewidmet. Aus- 
gehend von den natürlichen Zahlen werden durch 
Symmetrisierung von Addition und Multiplikation die 
ganzen und die rationalen Zahlen eingeführt und an- 
schließend die wichtigsten algebraischen Strukturen 
(Gruppen, Ringe, Körper, Vektorräume) definiert. 


Ed Ze 
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Vektorräume und lineare Abbildungen zweier Vektor- 
räume aufeinander stehen im Mittelpunkt der Be- 
trachtungen dieses Teiles. 

Die grundlegenden Begriffe in der Darstellung der 
Aralysis im dritten Teil sind Cauc#vfolge und gleich- 
mäßige Konvergenz. Die Eigenschaften der aufeinem 
Intervall der reellen Achse definierten stetigen reellen 
Funktionen und der Regelfunktionen (das sind die 
gleichmäßigen Limites von Treppenfunktionen) wer- 
den ausführlich untersucht und an Hand dieser Bei- 
spiele der Begriff des normierten vollständigen Vektor- 
raumes eingeführt. Modern ist auch die Darstellung 
der Integrationstheorie: Das Integral wird zunächst 
definiert auf dem Vektorraum der Treppenfunktionen 
über einem kompakten Intervall und anschließend 
durch stetige Fortsetzung erweitert auf den Raum der 
Regelfunktionen. Der R” ist konsequent als metri- 
scher Raum aufgefaßt. So ergeben sich die Aussagen 
über Vektorfunktionen bzw. Funktionen mehrerer 
reeller Veränderlicher als natürliche Verallgemeine- 
rungen der entsprechenden Sätze bei reellen Funk- 
tionen einer reellen Veränderlichen. Am Ende des 
dritten Teiles findet sich eine kurze Einführung in die 
Theorie der Linienintegrale und der mehrfachen Inte- 
grale. Man vermißt hier jedoch die Definition des 
Inhaltes und die Berechnung von ‘ gekrümmten 
Flächen. 

Im vierten Teil (‚Mathematische Hilfsmittel und 
Methoden“) werden unendliche Reihen, Entwicklun- 
gen von Funktionen, die Integration auf nicht kom- 
pakten Intervallen sowie die einfachsten Typen ge- 
wöhnlicher Differentialgleichungen behandelt. 

Die Vorzüge des Buches liegen in seiner Gründ- 
lichkeit und der Klarheit der Darstellung. }Leider 
enthält es keine Übungsaufgaben. Trotzdem bietet 
es dem Studenten nach Meinung des Referenten eine 
sehr gute Vorbereitung auf tiefergreifende mathe- 
matische Studien. 


Dresden H. LANGER 


Jahrbuch 1958 der Wissenschaftlichen Ge- 
sellschaft für Luftfahrt (WGL). Herausg. von 
H. Blenk. 300 S. m. zahlr. Abb. Braunschweig 1959. 


. Friedr. Vieweg & Sohn. Preis geb. 38,— DM. 


Das in jedem Jahr erscheinende WGL-Jahrbuch 
ist nicht nur für die Teilnehmer an der jeweiligen 
Tagung eine willkommene Zusammenstellung der dort 
gehaltenen Vorträge (einschl. Diskussionen), sondern 
gibt auch einem großen Kreis von interessierten Wis- 
senschaftlern und Ingenieuren, die an der Tagung 
nicht teilgenommen haben, einen aufschlußreichen 
Überblick über den neuesten Stand der Luftfahrtfor- 
schung und die gegenwärtig im Vordergrund stehen- 
den Forschungsthemen. Der vorliegende Band bringt 
von den meisten der 45 auf der Stuttgarter Tagung 
(8.—11. 10.1958) gehaltenen Vorträgen den Wort- 
laut, von den übrigen wenigstens eine ausführliche 
Übersicht. Der Festvortrag von S. Rurr galt der 
Frage ‚Flugzeug und Mensch“. Er wurde ergänzt 
durch einen Vortrag „Über die Einheit von Aus- 
rüstung und Flugzeug‘ von H. Kopp. Den Leser der 
ZAMM interessieren vor allem die Vorträge über Un- 
tersuchungen, die mit der Tragflügeltheorie zu- 
sammenhängen (K. Gersten, Tragflügeltheorie bei 
Unterschallgeschwindigkeit; F. KowALKE, Die flug- 
mechanischen Beiwerte von Tragflügeln bei Unter- 
schallgeschwindigkeit; F. WEGENER, Zur Program- 
mierung von Verfahren für die Berechnung der Auf- 
triebsverteilung an Tragflügeln; F. TmoMAs, Aero- 
dynamische Eigenschaften von Pfeil- und Deltaflügeln 
in Bodennähe), über Fragen der Gasdynamik 
(K. OswArıtscH und I. TEıreL, Die Pulsationen von 
Stoffdiffusoren ; M. SCHÄFER, Zur Stabilität von stoß- 


freien transsonischen Gasströmungen; A. NAUMANN 
und J. Prausz, Strömungsvorgänge im Rohreinlauf 
bei hohen Geschwindigkeiten), über Grenzschicht- 
erscheinungen und Konvektionsströmungen 
(E. Becker, Instationäre Grenzschichten hinter Ver- 
diehtungsstößen und Expansionswellen; W. PEcHAU, 
Ein Näherungsverfahren zur Berechnung der ebenen 
und rotationssymmetrischen turbulenten Grenzschicht 
mit beliebiger Absaugung oder Ausblasung; N. ScHoLz, 
Zur rationellen Berechnung laminarer un turbulenter 
kompressibler Grenzschichten mit Wärmeübertra- 
gung), über Gitterströmungen (K.H. GrEwE, 
Druckverteilungsmessungen an ebenen Schaufel- 
gittern bei hohen Unterschallgeschwindigkeiten; 
J. REHBAcH, Experimentelle Untersuchung der Spalt- 
verluste in Turbinenleiträdern mit bewegter Spalt- 
wand; G. MUESMANN, Zusammenhang der Strömungs- 
eigenschaften des Laufrades eines Axialgebläses mit 
denen eines Einzelflügels; W.-H. Isay, Die Strömung 
durch ein ‚schwingendes und rotierendes radiales 
Schaufelgitter), über Schwingungen in Triebwerken 
(G. HırscH, Beitrag zur Untersuchung der Biege- 
schwingungen von umlaufenden Triebwerk-Lauf- 
schaufeln mit nachgiebiger Einspannung) und über 
Fragen der Regelung (R. WoTscHkE, Grundsätz- 
liche Betrachtungen zur Flugregelung; E. MEYER- 
HarrwıG, Zwei neue elektronische Bauelemente für 
die Regeltechnik, der Linearkompensator und der 
Drahtheißleiter). In das Gebiet der Unternehmensfor- 
schung gehört außer dem Vortrag W. SPILGER, The 
application of weapon systems analysis and military 
operations research techniques to military planning, 
im Grunde auch der Vortrag von G. BAoKHAUS: 
Grundbeziehungen für den Entwurf optimaler Ver- 
kehrsflugzeuge. Vorträge über Meßverfahren und 
Meßtechnik, Konstruktionsfragen, über Werkstoffe 
usw. vervollständigen das reichhaltige und farbige 
Bild, das auch dieser (wie immer vorzüglich ausge- 
stattete) Band des WGL-Jahrbuches von der Luft- 
fahrtforschung entwirft. 


Dresden H. HEINRICH 


W.Matz, Wärme- und Stoffaustausch in 
Arbeitsdiagrammen auf projektiver Grund- 
lage. VIII + 142 S. m. 44 Aufg. u. 56 Abb. Berlin/ 
Göttingen/Heidelberg 1960. Springer-Verlag. Preis 
geb. DM 24,—. 


In dem vorliegenden Buch hat der Verfasser erfolg- 
reich den Versuch unternommen, das vollständige 
Viereck der Geometrie der Lage mit seinen Beziehun- 
gen zu vier harmonischen Punkten auf den Wärme- 
und Stoffaustausch in verfahrenstechnischen Pro- 
zessen anzuwenden. Nach einer kurzgefaßten Ein- 
führung in den analytischen und konstruktiven Zu- 
sammenhang zwischen Viereck und harmonischen 
Punkten und in die Analogie zwischen Kräften und 
Kraftrichtungen einerseits und strömenden Massen 
und deren Strömrichtungen andererseits werden die 
in Wärmetauschern, bei der Rektifizierung und bei 
der Flüssig-Flüssig-Extraktion vorliegenden Aus- 
tauschprobleme geometrisch ausführlich behandelt. 
Es ist erstaunlich und für den Mathematiker wie für 
den Verfahrensingenieur gleichermaßen interessant, 
mit welcher Klarheit, Anschaulichkeit und verblüffen- 
den Einfachheit die physikalischen Vorgänge sich 
durch diese erstmalig verwendete grafisch-geome- 
trische Methode zur Darstellung bringen lassen. 

Der Stoff ist ohne weitergreifenden analytischen 
Aufwand logisch entwickelt und klar und leicht faßlich 
dargeboten, vertieft und belebt durch 33 der prak- 
tischen Verfahrenstechnik entnommene Beispiels- 
aufgaben. Das Buch verdient stärkste Beachtung 
sowohl seitens der Mathematiker als treffliches Schul- 
beispiel als auch seitens der Studierenden und prak- 
tisch tätigen Ingenieure der Verfahrenstechnik. 


Dresden H. FALtın 


Festkörpermechanik zu 1 h, 

Der erste Band der geplanten Buchreihe lie 
mehr vor; er enthält ausschließlich theoretische i 
sätze. Das erste Rapid, verfaßt von 8. C. Hunter, 
behandelt das Problem der Wellenausb: 
visco-elastischen Körper. Nach einer phy 
Deutung des FIRTNERGSDEN, Verhaltens werc 
Gleichungen für Spannung, ıngsgesc 
keit, Dehnung und Dehnungsgeschwindigkeit 
siert, so daß die Anwendung der Lartacr-Tı j 
mation, der Foukrer-Transformation oder einer Ope- 
ratoren-Methode möglich wird. Anschließend wird 
eine Gruppe von Lösungen untersucht, welche für die 
Ausbreitung von Störungen maßgebend sind. Die 
gefundenen Ergebnisse werden Versuchsergebnissen 
gegenübergestellt. Eine tensorielle Darstellung der 
Ausgangsgleichungen im dreidimensionalen Falle, 
sowie Bemerkungen zur Anwendbarkeit der dyna- 
mischen Visco-Elastizität auf Probleme des Ultra- 
schalls ergänzen die Abhandlung. 


Der zweite Artikel, verfaßt von K. MARGUERRE, 
behandelt die Anwendung der Matrizenrechnung auf 
Probleme der Stabschwingungen. Nach Aufstellung 
der Grundgleichungen werden die einzelnen Glieder 
der auftretenden Matrizen für einen Spezialfall ange- 
geben. Ferner werden Probleme mit stückweise kon- 
stanter Masse pro Längeneinheit vorgeführt. Weitere 
Abschnitte behandeln Besonderheiten, wie sie durch 
Einspannungen, Auflagerbedingungen sowie Einzel- 
massen hervorgerufen werden. 

Die dritte Arbeit, verfaßt von H.G. Horkıns, 
betrifft die dynamische Ausbildung von Hohlräumen 
in Festkörpern, besonders in zähen Metallen. Es 
liegen sowohl rein elastische als auch elastisch- 
plastische Untersuchungen vor, wobei die letzteren 
heute von größerem Interesse sind. Eine außerordent- 
lich große Mannigfaltigkeit von Sonderproblemen in 
Physik und Technik hängen mit diesem Problem eng 
zusammen. Ähnliche Probleme treten z.B. bei der 
Belastung von Behältern durch Innendruck auf; desgl. 
bei der Autofrettage von Geschützrohren. Der Verf. 
beschränkt sich auf die Untersuchung unter folgenden 
Annahmen: Die Bewegung ist kugelsymmetrisch; die 
Bewegung der Gasblase ist adiabatisch. Das Metall 
verhält sich inkompressibel und befolgt im elastischen 
Gebiet das Hookzsche Gesetz, im plastischen Gebiet | 
die TrescAasche Fließbewegung. Verfestigungs- und 
Dehnungsgeschwindigkeitseinflüsse werden empirisch 
berücksichtigt. Die theoretischen Grundgleichungen 
werden für diese Voraussetzungen diskutiert, insbe- 
sondere auch für große Deformationen und mit Be- 
rücksichtigung der bei einer plötzlichen Schockein- 
wirkung (Explosion) auftretenden Trägheitsglieder. 


Ein Aufsatz von W. T. Korter behandelt die Grund- 
lagen der elastisch-plastischen Theorie. Die Spezial- 
fälle des vollkommen plastischen Materials und des 
Materials mit Verfestigung werden für kleine Ver- 
formungen diskutiert. Das Eindeutigkeitsproblem 
wird mit Berücksichtigung der Spannungsgeschwin- 
digkeit nachgewiesen. Die Minimalprinzipien werden 
für Spannungen und Dehnungen, sowie für Spannungs- 
und Dehnungsgeschwindigkeit abgeleitet. Weitere 
Abschnitte behandeln die Grenzlastanalyse und den 
plastischen Bruch, sowie einige Existenzbeweise. 
Ein Abschnitt von W.A. Grren behandelt die 
Dispersion bei elastischen Wellenbewegungen in Stä- 
ben. Nach Erörterung der klassischen exakten Lö- 
sungen werden die bekannten Näherungsmethoden 
beschrieben und einander gegenübergestellt. Schließ- 
lich werden die Dispersionserscheinungen bei der Aus- 


- 
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 breitung hochfrequenter Wellen in Stäben am Beispiel 
des Kreiszylinders diskutiert. 


Im sechsten Abschnitt behandelt P. CHapwick 
dynamische Probleme der Wärmespannungen. Nach 
Aufstellung der Grundgleichungen für nicht umkehr- 
bare thermodynamische Vorgänge im elastischen 
Festkörper werden zunächst die ebenen harmonischen 
thermoelastischen Vorgänge untersucht und mit Ver- 
suchsergebnissen verglichen. Die Erörterung von 
Randwertproblemen erfolgt am Beispiel des Kreis- 
zylinders, wobei thermoelastische RavLeıcHwellen 
und Schwingungen berechnet werden. Weitere Über- 
legungen beziehen sich auf kugelsymmetrische und 
achssymmetrische Wellenbewegungen. 


B. A. Bıtgy behandelt stetige Verteilungen von 
Dislokationen. BURGERSvektor und Torsionstensor 
führen zu einer Modifikation das Spannungstensors 
gegenüber der bisherigen Kontinuumsmechanik, wel- 
che mit einer Veränderung des Kristallgitters in Zu- 
sammenhang gebracht wird; auf diese Weise lassen 
sich Einflüsse bleibender Deformationen und Eigen- 
spannungen berücksichtigen. 


Ein letztes Kapitel von R. MukI beschäftigt sich 
mit Problemen des elastischen Halbraumes und der 
dieken Platte. Es wird die Wirkung einer Kraft- 
einleitung, herrührend von einem starren Körper, 
untersucht. Ein weiteres Problem behandelt die 
Wärmespannungen in einem halbunendlichen Körper 
undin einer dicken Platte bei stationärer Temperatur- 
verteilung. 


München H. NEUBER 


S. Karlin, Mathematical Methods and Theory 
in Games, Programming, and Economics. 
Vol. I-+II. X + 4338. und XI + 386 S. London/ 
Paris 1959. Pergamon Press. Preis geb. pro Band 
75 s. net. 


Der Titel des vorliegenden Werkes charakterisiert 
insbesondere den ersten Band ganz ausgezeichnet; 
Die ersten 4 Kapitel dieses Bandes sind der Spiel- 
theorie, die nächsten 3 den linearen Programmen und 
die letzten 2 Kapitel der mathematischen National- 
ökonomie gewidmet. Die Darstellung ist nicht auf 
größtmögliche Allgemeinheit ausgerichtet, aber sehr 
ausgefeilt und durchdacht. Konzessionen an die 
mathematische Strenge werden nicht gemacht. Viel- 
fältig ist neues und originelles Material in das Werk 
eingearbeitet worden. Bei der Behandlung der Spiel- 
theorie werden nur Zweipersonenspiele betrachtet. 
Der Stoff wird meist in einer Form dargeboten, welche 
unmittelbare Schlüsse auf Verallgemeinerungen zu- 
läßt. So wird für den Minimaxsatz für endliche Zwei- 
personenspiele auch ein Beweis mittels des Fixpunkt- 
satzes von KAKUTANI gegeben, so daß also im Prinzip 
der Satz für lineare topologische Räume bewiesen ist. 
Es ist interessant zu sehen, daß die Kapitel über Spiel- 
theorie im ersten Band mit den entsprechenden Ka- 
piteln des beinahe zur selben Zeit erschienenen deutsch- 
sprachigen Buches von BURGER, das methodisch und 
inhaltlich — allerdings in viel engerem Rahmen — 
eine ganz ähnliche Zielsetzung hat, sehr viel gemein- 
sam haben. Die Kapitel über die Programme geben 
neben der Theorie der linearen Programme auch die 
wichtigsten Rechenmethoden. Besonderes Interesse 
verdient natürlich die hier gegebene und, wie mir 
scheint, erste lehrbuchmäßige Darstellung nicht- 
linearer, insbesondere konvexer Programme. Als 
wesentliches mathematisches Werkzeug wird eine von 
FEncHEL entwickelte Theorie konvexer Funktionen 
benützt, die vom bekannten Begriff der Polarität aus- 
geht. Spieltheorie und mathematische Programme 
haben stark das Gesicht der modernen mathemati- 
schen Nationalökonomie beeinflußt. Es erscheint da- 
her methodisch zweckmäßig, dieses Gebiet nach der 
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Spieltheorie und den Programmen zu behandeln. 
Zunächst werden verschiedene Aspekte der Gleich- 
gewichtstheorie betrachtet, die schließlich im Modell 
von ARROW-DEBREU gipfeln. Dieses Kapitel enthält 
auch eine gut lesbare Einführung in die Theorie der 
positiven Matrizen. Das letzte Kapitel des ersten 
Bandes beschäftigt sich mit der Welfare-Economics, 
die jaim Gegensatz zur klassischen Nationalökonomie 
fast ausschließlich von angelsächsischen Gelehrten 
entwickelt wurde, und enthält auch das klassische 
von NEUMANNsche Modell einer expandierenden Wirt- 
schaft und dessen Verallgemeinerungen. Jedes Ka- 
pitel enthält am Schlusse Hinweise historischer Art 
und Verweise auf das außerordentlich reichhaltige 
Literaturverzeichnis am Schluß des Bandes. Der An- 
sicht des Autors: “his (WALRASs) theory of general 
equilibrium did for economics what Newton did for 
physies“ wird man allerdings nur schwer zustimmen 
können. Jedem Kapitel ist auch eine Reihe von Pro- 


 blemen beigegeben, wie es ja heute bei amerikanischen 


Büchern fast allgemein üblich ist. Drei Anhänge be- 
schließen den ersten Band. Der erste enthält eine 
Einführung in die lineare Algebra, der zweite beschäf- 
tigt sich mit konvexen Mengen und Funktionen, und 
der dritte enthält u.a. Fixpunktsätze und den Be- 
griff der Wahrscheinlichkeitsverteilung. Während 
der erste Band ausschließlich endliche Zweipersonen- 
spiele behandelt, beschäftigt sich der zweite Band mit 
unendlichen Zweipersonenspielen. Er ist zwar als 
Fortsetzung des ersten Bandes gedacht, aber doch 
völlig unabhängig von diesem lesbar und stellt ein in 
sich abgeschlossenes Werk dar. Bein äußerlich kommt 
das auch dadurch zum Ausdruck, daß das erste ein- 
leitende Kapitel des ersten Bandes über die Definition 
eines Spieles und den Begriff der Minimax-Theoreme 
unverändert im zweiten Band nochmals abgedruckt 
ist. Der zweite Band enthält auch dieselben Anhänge 
und dasselbe Literaturverzeichnis. Dieser Band ist 
weniger systematisch angelegt als der erste, aber das 
liegt in der Natur der Sache. Der Autor, der sich 
selbst sehr aktiv mit unendlichen Spielen beschäftigt 
hat, gibt eine interessante Begründung dafür: “sol- 
ving games is like solving differential equations, a 
combination of good guesswork ...“. Der etwas fort- 
geschrittene Leser wird mit Interesse zu diesem Band 
greifen. Zunächst werden allgemeine Klassen von 
Spielen behandelt: Spiele, deren Auszahlungsfunktionen 
Polynome sind, und dann solche, deren Auszahlungs- 
funktionen konvex sind. Die weiteren 5 Kapitel sind 
einer Reihe von interessanten Beispielen gewidmet. 
Wir erwähnen insbesondere das Kapitel über die Ver- 
allgemeinerung des NEYMAN-PEARSoNschen Lemmas, 
dessen Schlußweite auf eine Reihe von Spielen ange- 
wendet wird, die nicht über dem Einheitsquadrat ge- 
spielt werden. 


Wien L. SCHMETTERER 


0. Tietjens, Strömungslehre. Physikalische 
Grundlagen vom Technischen Standpunkt. Band]. 
Hydro- und Aerostatik. Bewegung der idealen 
Flüssigkeit. XVI + 536 S. m. 496 Abb. Berlin/Göt- 
tingen/Heidelberg 1960. Springer-Verlag. Preis geb. 
66,— DM. 


Der Name Tiersess ist bekanntgeworden durch 
das erstmalig 1929 im gleichen Verlag erschienene 
Buch Hydro- und Aeromechanik, womit der Verfasser 
die verdienstvolle Aufgabe übernommen hatte, nach 
L. PRANDTLs Göttinger Vorlesungen ‚„Hydrodynamik 
und Aerodynamik‘ sowie ‚„‚Aeromechanik‘ ein zwei- 
bändiges Lehrbuch zu schreiben. Neben der vorlie- 
genden neuen Bearbeitung wirkt der 238 Seiten um- 
fassende erste Band des alten Lehrbuches wie eine 
stark gekürzte Zusammenfassung des jetzt neu ge- 
botenen. Die sehr breiten, gemächlichen, ja gemüt- 
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lichen Beschreibungen atmen Ruhe. Man kann hier 
von einer konservativen Darstellung der Strömungs- 
lehre sprechen. Durch viele Bilder belebt, führt der 
Verfasser den Leser in die Grundlagen ein. Er spricht 
über den Flüssigkeitsdruck, die Stabilität von Schiffen 
und Ballons, Probleme der Erdatmosphäre, Kapil- 
larität, über ebene und drehsymmetrische Potential- 
strömungen, Flüssigkeitswellen und Tragflügeltheorie. 
Die Probleme der Gasdynamik und Reibung sind 
einem geplanten zweiten Band vorbehalten. 


Die Grundbegriffe Zähigkeit, Dichte und derglei- 
chen werden sehr ausführlich erläutert, jedoch könnte 
hier vieles beim Leser als schon von der Oberschule 
her bekannt vorausgesetzt werden. 


Empfehlenswert wäre auch, daß das Buch in kom- 
menden Auflagen das reformierte Maßsystem benutzt 
nach den Empfehlungen der 10. Generalkonferenz für 
Maß und Gewicht aus dem Jahre 1954, dem sich viele 
Autoren anderer Lehrbücher dieses Fachgebietes 
schon angeschlossen haben, wobei besonders wichtig 
ist, daß als Krafteinheit, wenn nicht gleich das New- 
ton, so doch das Kilopond benutzt wird und das 
Kilogramm für die Masseneinheit reserviert bleibt. 


Als besonders nützlich für den Lernenden muß die 
exakte Darstellung von Strom- und Äquipotential- 
linien vieler Strömungsfelder bezeichnet werden. Die 
Ablesung von Koordinaten einzelner Punkte aus 
diesen Feldern wird erleichtert durch symmetrische 
Skalen, die rechts und links oben und unten an den 
Zeichnungen angebracht sind, und eine mitgelieferte 
durchsichtige Scheibe, auf die ein Liniennetz gezeich- 
net ist. 

Bei den Anwendungen hat weniger die Suche nach 
aktuellen Beispielen im Vordergrund gestanden, als 
die Absicht, historisch bedeutungsvolle und gesicherte 
Ergebnisse zu vermitteln. So wird auf 18 Seiten der 
Leser in die gasdynamischen und meteorologischen 
Probleme des Freiballons eingeführt. Ebenfalls wird 
die FUHRMANNsche Dissertation mit dem Aufbau von 
Drehkörperformen aus Quellen und Senken, die für 
die Gestaltung von Luftschiffkörpern von Bedeutung 
war, eingehend referiert. 

Die Reproduktionen von Strömungsfotografien 
ähnlich denen aus bekannten Göttinger Strömungs- 
filmen sind nützlich, weil gerade sie von hervorragen- 
der pädagogischer Bedeutung für den Lernenden sind. 


Eu a a a ih 


Es werden u.a. dargestellt die Bildung des Anfahr- 
wirbels mit 10 Fotos verschiedener Entwicklungs- 
phasen. “ 

Es ist ein Buch von besonderer Eigenart, auf dessen 
Anschaffung weder der Lernende noch der Fachmann 
verzichten sollten. j 


Dresden W. ALBRING 


K. Klöppel und J. Scheer, Beulwerte ausgesteif- 
ter Rechteckplatten. Kurventafeln zum direkten 
Nachweis der Beulsicherheit für verschiedene Steifen- 
anordnungen und Belastungen. 107 S. m. 45 Abb. u. 
103 Beulwerttafeln. Berlin 1960. Verlag von Wilhelm 
Ernst & Sohn. Preis geb. DM 44,—. 


Die Forderungen des Leichtbaus haben zwangs- 
läufig zur Verwendung dünner ausgesteifter Bleche 
und Platten geführt, für deren statische Berechnung 
ein großes wissenschaftliches Schrifttum entstanden 
ist. Auch die zuständigen Normungsvorschriften 
(DIN 4114) lassen erkennen, wie umfangreich und wie 
kompliziert die vom Statiker vorzunehmenden Be- 
rechnungen in der Regel sein müssen, um allen Bedin- 
gungen gerecht zu werden. Um diese große Arbeit 
des Statikers abzukürzen, haben die Verfasser sich 
die Aufgabe gestellt, ein Tafelwerk herauszugeben, 
das dem-direkten Nachweis der Beulsicherheit dient, 
wobei verschiedene Belastungsfälle und verschiedene 
Anordnungen der Versteifungen berücksichtigt sind. 
Die umfangreiche Arbeit wurde mit Unterstützung 
der Deutschen Forschungsgemeinschaft im Institut 
für Praktische Mathematik der Technischen Hoch- 
schule Darmstadt unter Hinzunahme eines programm- 
gesteuerten Digitalrechners (IBM 650) durchgeführt 
(später stand auch der schnellere IBM 704 zur Ver- 
fügung). Das Buch enthält ausgezeichnete Kurven- 
tafeln mit klarer Angabe der jeweils erforderlichen 
Erläuterungen, so daß der gesamte Leichtbau, ins- 
besondere Stahlbrückenbau, Stahlhochbau und Kran- 
bau das Werk zur Anschaffung für die statischen 
Büros dankbar begrüßen werden. Ebenso wird das 
Buch an allen einschlägigen Lehrstühlen der Tech- 
nischen Hochschulen sowie Forschungsstellen der 
Industrie starkes Interesse finden. 


München H. NEUBER 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


A. Linder, Statistische Methoden für Natur- 
wissenschaftler, Mediziner und Ingenieure. Dritte 
Aufl. 4848. Basel und Stuttgart 1960. Birkhäuser 
Verlag. Preis geb. 58,60 DM. 


R. Sauer, Einführungin dietheoretische Gas- 
dynamik. Dritte verbesserte Auflage. XII -+ 2148. 
m. 120 Abb. Berlin/Göttingen/Heidelberg 1960. 
Springer-Verlag. Preis geb. DM 29,70. 


Ya.L. Geronimus, Polynomials Orthogonal 
on a Cirele and Interval. IX + 210 S. m. 9 Ta- 
feln. Oxford/London/New York/Paris 1960. Perga- 
mon Press. Preis geb. 60 s.net. 


G.B. Thomas, Jr. (Department of Mathematics, 
Massachusetts Institute of Technology), Caleulus 
and Analytie Geometry. Third Edition. XII + 
1010 S. London 1960. Addison-Wesley Publishing 
Company, Inc. Preis geb. $ 9,75. 


W.Heyn,Stichprobenverfahren in der Markt- 
forschung. 108 8. und 7 Tafeln. Würzburg 1960. 
Physica-Verlag. Preis geb. DM 22,—. 


H. Schlitt, Systemtheorie für regellose Vor- 
gänge. Statistische Verfahren für die Nachrichten- 
und Regelungstechnik . XIT + 344 S. m. 167 Abb. 
Berlin/Göttingen/Heidelberg 1960. Springer-Verlag. 
Preis geb. DM 49,50. 


K. Jacobs, Neuere Methoden und Ergebnisse 
der Ergodentheorie. (Ergebnisse der Mathematik 
und ihrer Grenzgebiete, Neue Folge, Heft 29). VI -+ 
214 S. Berlin/Göttingen/Heidelberg 1960. Springer- 
Verlag. Preis brosch. DM 49,80. 


J. Palacios, Analyse Dimensionnelle. Traduit 
de P’espagnol par J. Prevot. XV + 206 8. Paris 1960. 
Gauthier-Villars. Preis brosch. $ 8,50. 


G. Birkhoff, Hydrodynamics. A Study in Logie, 
Fact, and Similitude. Second edition, revised and en- 
larged. 184 S. m. 28 Abb. Princeton 1960. Princeton 
University Press. Preis geb. $ 6,50. 


W. Späth, Zahl — Maß — Bild. Grundfragen der 
Meßtechnik. 227 S. m. 129 Abb. Stuttgart 1960. 
Verlag A. W. Gentner KG. Preis geb. DM 29,50. 
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F. Bosnjakovid, Technische Thermodynamik. 
II. Teil. 3. Auflage. (Wärmelehre und Wärmewirt- 


' schaft in Einzeldarstellungen, Band 12). XX + 
480 S. m. 359 Abb. u. Aufgabensammlung. Dresden 


und Leipzig 1960. Verlag von Theodor Steinkopff. 
Preis geb. DM 23,40. 


Kybernetik. Zeitschrift für Nachrichtenüber- 
tragung, Nachrichtenverarbeitung, Steuerung und 
Regelung im Organismus und in Automaten. Band I, 
Heft 1. 56 8. m. 72 Abb. Berlin/Göttingen/Heidel- 
berg 1961. Springer-Verlag. Preis DM 12,80. 


L. M. Milne-Thomson, Plane Elastic Systems. 
VIII + 221 S. m. 76 Abb. Berlin/Göttingen/Heidel- 
berg 1960. Springer-Verlag. Preis brosch. DM 49,80. 


H. P. Künzi, Quasikonforme Abbildungen. 
(Ergebnisse der Mathematik und ihrer Grenzgebiete, 
Neue Folge, Heft 26.) VIII + 182S.m. 35 Abb. Berlin/ 
Göttingen/Heidelberg 1960. Springer-Verlag. Preis 
brosch. DM 39,—. 


K.-R. Biermann, Vorschläge zur Wahl von 
Mathematikern in die Berliner Akademie. 
Ein Beitrag zur Gelehrten- und Mathematikerge- 
schichte des 19. Jahrhunderts. (Abhandlungen der 
Deutschen Akademie der Wissenschaften zu Berlin, 
Klasse für Mathematik, Physik und Technik, Jg. 
1960, Nr. 3). 75 8. Berlin 1960. Akademie-Verlag. 
Preis brosch. DM 7,50. 


6. Köthe, Topologische lineare Räume I. 
XII + 456 S. Berlin/Göttingen/Heidelberg 1960. 
Springer-Verlag. Preis geb. DM 78,—. 


R. Bellman, Adaptive Control Processes. 
XVI + 255 S. Princeton, N. J., 1961. Princeton 
University Press. Preis geb. $ 6.50. 


A. Kratzer und W. Franz, Transzendente 
Funktionen. (Mathematik und ihre Anwendungen 
in Physik und Technik, Reihe A, Band 28.) XIII + 
375 8. m. 58 Abb. u. 1 Tab. Leipzig 1960. Akademi- 
sche Verlagsgesellschaft Geest & Portig K.-G. Preis 
geb. DM 39,—. ö 


A.Duschek, Vorlesungen über höhere Mathe- 
matik. Zweiter Band: Integration und Differen- 
tiation der Funktionen von mehreren Veränderlichen. 
Lineare Algebra. Tensorfelder. Differentialgeometrie. 
Zweite, neu bearbeitete Aufl. VIII + 401 S. m. 
a Wien 1958. Springer-Verlag. Preis geb. 


K. Löffler, Die Berechnung von rotierenden 


Scheiben und Schalen. VIII + 241 S. m. 92 Abb., 
4 Kurvenblättern und 4 Formblättern als Anlage. 
Berlin/Göttingen/Heidelberg 1960. Springer-Verlag. 
Preis geb. DM 42,—. 


E.Kähler, Innerer und äußerer Differential- 
kalkül. (Abhandlungen der Deutschen Akademie der 
Wissenschaften zu Berlin, Klasse für Mathematik, 
Physik und Technik, Jg. 1960, Nr. 4) 32 S. Berlin 
1960. Akademie-Verlag. Preis brosch. DM 5,50. 


V. Fock, Theorie von Raum, Zeit und 
Gravitation. XIX + 501 8. Berlin 1960. Akademie- 
Verlag. Preis geb. DM 42,50. 


NACHRICHTEN 


The Sixth Congress on Theoretical and Applied Mechanics 
and the Symposium on 
“High Speed Computation Methods and Machines” 
Delhi, India 


The Indian Society of Theoretical and Applied Me- 
chanies held their Sixth Congress in Dehli this year 
between 24. 12. 60 and 26. 12. 60 under the auspices 
of the University of Dehli. It was preceded by a one- 
day symposium on High Speed Computation Methods 
and Machines arranged, by the Society jointly with 
the International Business Machines Corporation. 
About 150 delegates and members, including those 
from Burma, France, Germany, Pakistan, Poland, 
U.K. and U. S$. A. attended the symposium and the 
Congress. 

After the welcome adress by the Vice-chancellor of 
the University of Dehli who spoke about the need 
for University institutions contributing fully to the 
development of a Society of this type in the impor- 
tant branches of Mechanics, Dr. D. S. KoTHARt in- 
augurated the symposium with a short address devo- 
ted to the interaction of micro and macromechanics 
on the one hand and the need for a greater application 
of computing and controlling devices in respect of 
space travel and communication on the other hand. 
Dr. A.N. KnostA, the President of the Congress 
reviewed the work of the Society during the past five 
years, particularly from the point of view of the con- 
tribution of various schools of research activity in 
the country. A. K. CHAUDHURI in his half-hour ad- 
dress gave a review of the current position in the 
development and application of high speed comput- 
ing machines and their limitations and emphasised 
the advances taking place in pure and applied mathe- 


matics as a result of the advent of high speed com- 
puting devices. 

R. SAUER (Germany) in his half-hour address on 
experimental mathematics on hish speed computers 
reported on the error analysis for difference methods 
and numerical investigation of the stability of shock 
free transonic flows past profiles. R. NARASIMHAN 
described the essential features of a digital computer 
being built and operated by the Atomic Energy Com- 
mission of India. M. K. JAın spoke on various collo- 
cation methods based on TCHEBYSCHEFF’sS norm to & 
wide class of problems including non-linear boundary 
value problems. G. Pıckerr (Wisconsin, U.S.A.) 
described the problems under investigation at the 
University of Wisconsin digital computer centre. 
B. M. BELGAUMKAR reported on application of RLC 
network models for the study of stresses in crank 
shafts. S.N.B. Murr#Y described some analogue 
computer studies on the stability of single and two 
shaft power generating systems, with particular re- 
ference to gas turbine power plants. J. K. SEN repor- 
ted on the use of an analogue computer in the analysis 
of hydraulic surge problems and compared the results 
with model tests. S. KrRısHNAn gave details of the 
high speed data handling systems for the flow down 
wind tunnel operating at Bangalore. H. B. MoHANTI 
spoke on the criterion of stability of feed back systems 
with statistical inputs. A. K. CHAUDHURI reported 
on the stability of a third order non-linear servo-me- 
chanism with field built up time and hysteresis. 


| 
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Sixth Congress 


The half-hour addresses were delivered by R. SAUER 
(Germany), J. Kampr DE Ferıer (France), H. ZoRSKI 
(Poland), H.B. Sauırn (England), H. LangHaaR 
(U.S.A.) and G. N. MESHOHERYAKOV (USSR.). 

R. SAvEr spoke on three dimensional problemsin gas 
dynamics and showed the application of characteri- 
sties grid methods to unsteady planar and steady spa- 
tial gas dynamic problems. J. KamrE DE FeRIET 
described the statistical mechanices of linearized plane 
sound waves in a perfect gas filling the whole space. 
H. ZoRsKI gave a review of Polish papers and some 
new results on plates with discontinuous boundary 
conditions. H.B. Saure spoke on the theory of 
radial jet flows and described experimental work con- 
cerned with such flows. H. LAn@HAAR gave an account 
of the general theory of stability of conservation sy- 
stem. G. N. MESHCHERYAKOV spoke on the nature of 
self-exeited vibration in metal cutting processes. 

Fluid Dynamics: R. BaALLABH spoke on the un- 
steady exact solution of the NAvIER-STOKES equa- 
tions. M.N.L. NAarAasımHAan described the laminar 
non-Newronian flow in a porous pipe and showed 
that the velocity profile and pressure gradient are 
considerably affected. M. K. JAIN, using second order 
stress strain velocity relations discussed the effect of 
cross-viscosity in fluid flow problems. S. N. SıycH 
and G. A. NARIBOLI reported their work on the flow 
of conducting fluids in circular pipes under transverse 
magnetic field. A. S. GuprTA spoke on the unsteady 
free convection flow of a compressible fluid past an 
accelerated flat plate subjected to suction. R. S. NAn- 
DA discussed the heat transfer by laminar flow be- 
tween torsionally oscillating disks. C. N. KAUuL spoke 
on the growth and decay of week discontinuities in 
magnetohydrodynamics. G. L. SAINI reported on the 
theory of shock waves in relativistic magnetofluid- 
dynamics. L.V.K.V. SAarmA discussed the shear 
flow of a conducting fluid in the presence of a magne- 
tie field. S. K. SHARMA and A. Ü. SRIVASTAVA spoke 
on the torsional oscillations of a plane in a visco-elastie 
liquid. G. N. SAarmA described the steady boundary 
layer with suction using KARMAN integral equations. 
J. R. Rao discussed the ’strength‘ of the basie field 
equations of compressible flows. Miss S. L. RATHNA 
and Miss G. K. RAJESWARI spoke on the superposa- 
bilityin non-NEwronian fluids with variable viscosity. 
and cross-viscosity coefficients. 8S.N.B. Murruy 
discussed the high intensity shock waves taking ac- 
count of ionisation and radiation effects. J. N. KAPrur 
and his coworkers gave an account of their recent 
work in magneto-hydrodynamics. P. C. Jaın discus- 
sed the theory of hydromagnetie turbulence flows. 
R. N. BHATTACHARYA spoke on the characteristics of 
accelerated ship waves and graphs of navigational 
interest. P.N. CHATTERJEE and his coworkers de- 
scribe the scale effects in the model study of a hy- 
draulic surge tank. P. O. BuuraAnt described the non- 
equilibrium flows in magneto-gas-dynamics, 

Elasticity and others: G. A. NArısour and V. B. 
NYAYADHISH spoke on the wave propagation in semi- 
infinite eylindrical red. P.D.S. Verma and 8.B. 
SINHA described the propagation of symmetrical dis- 
turbancein a hypoelastie infinite plate of grade zero. 
D. N. MrrrA spoke on the indentation of a rigid para- 
boloid of revolution or a cone on a semi-infinite trans- 
versely isotropic medium. H.S. Paut and G. Parıa 
reported on magneto-elastic stresses in an infinite me- 
dium with a long eylindrical hole. J. Ramakanın 
discussed the finite strains in cireular eylinder and 
spherical shell under the action of radial body-forces. 


Herausgeber und Hauptschriftleiter: Prof. Dr.-Ing. habil. H. H einrich, Dresden A 27 
GmbH, Berlin W 8, Leipziger Straße 3—4; Fernsprecher: 22 0441. Telex- 
Heftes: 1009/41/5. Die Zeitschrift für Angewandte Mathematik und Mechanik 
Zuzüglich Bestellgeld. Doppelheft DM 10,—. Abbestellungen können nur bis 


tion the multiple valued solution for the k go 
a rectangular plate into right eylindrical can be 


ved beam. P. Purı and Avrar Sına# describe 
linear dynamie thermoelastic problem with a moving 
interface. 4 

E. Laxsumınarayana and B. Krsava Rao had 
solved the problem of bending of an aelotropie eircular 
plate into a spherical shell. U. A. Sasıry spoke on 
torsion of a circular eylinder reinforced by a square 
rod. P. N. Smrvkumar described the dislocation in 
material whose cross-section is bounded by two con- 
centric ellipses. 

K.G. ee P. R. Duar and K.N. Anar- 
THARAMAN discussed the hardness of low carbon steel 
at elevated temperatures. B.M. BELGAUMKAR spoke 
on the mechanics of singing propellers. J. N. Kapur 
and B. S. Jar discussed the modifications to ordinary 
multitubular charges to obtain better burning surfaces. 

A.K. Gayen and G. C. Roy spoke on the problem 
in multiple timeseries. A. BHATTAcHARYYA and R.K. 
GHosH discussed the analysis of chip formation me- 
chanism of restrieted contaet eutting tools. C.K. 
SARKAR described the effects of inclined eut-offs in 
reducing uplift pressures on the bases of hydraulie 
structures. 

Dr. D. S. Koruarı, Seientifie Adviser to the Mini- 
stry of Defence, has been elected the President and 
Prof. Ram BEHARI, (Dehli University) the Vice-presi- 
dent for the period 1961—1963. Prof. B.R. SerH 
(Kharagpur) continues as Secretary-Treasurer. The 
next congress in December 1961 will be held at the 
Indian Institute of Technology, Bombay. 

B. R. SETH 


Modern Developments in Heat Transfer 


Auf Wunsch des Veranstalters wird mitgeteilt, 
daß vom 14. bis 22. August 1961 an der Universität 
von Minnesota ein spezieller Sommerkursus über das 
im Titel genannte Thema unter der Anleitung von 
10 Wissenschaftlern aus den USA, England und 
Deutschland stattfindet. Die Teilnehmergebühr be- 
trägt $ 290 einschl. Unterbringung und Verpflegung. 
Interessenten wenden sich, bitte, an: Dr. E.R.G. 
ECKERT, Professor and Director, Heat Transfer 
Laboratory, University of Minnesota, Minneapolis 14, 
USA. 

In diesem Zusammenhang sei noch auf die Tagung 
„International Discussion on Heat Transfer“ auf- 
merksam gemacht, die vom 28. August bis 1. Septem- 
ber 1961 in Boulder stattfindet. 


IUTAM-Symposium Haifa 1962 


Im Frühjahr 1962 veranstaltet die IUTAM (Inter- 
national Union of Theoretical and Applied Mechanics) 
in Haifa/Israel ein Internationales Symposium unter 
dem Titel „Second Order Effects in Plasti- 
eity, Elastieity, and Fluid Dynamics“. Inter- 
essenten wenden sich wegen Einzelheiten bitte an 
Herrn Prof. Dr. H. GöRTLER, Freiburg/Br., Sonn- 
halde 90. Es werde aber schon jetzt darauf aufmerk- 
sam gemacht, daß Vortragsauszüge bereits bis zum 
15.77.61, die vollständigen Manuskripte bis zum 
1. 11. 61 einzureichen sind. 


27, Friedrich-Hegel-Str. 31. Verlag: Akademie-Verlag 
Nr. 011 773. Postscheckkonto: Berlin 350 21. Bestellnummer dieses 
erscheint monatlich. Bezugspreis: vierteljährlich DM 15,—. 


3 A ellu 4 Wochen vor Quartalsende anerkannt werden, andernfalls 
wird das folgende Quartal noch geliefert. Veröffentlicht unter Lizenznummer ZLN 5011 des Ministeriums für Kultur, 
Gesamtherstellung: VEB Druckerei ‚Thomas Müntzer‘“ Bad Langensalza (V/12/6) (1). Printed in Germany. 
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GE RHARD HEBER 


Vorlesungen 
über ausgewählte Kapitel der Quantenfeldtheorie 


(Abhandlungen der Deutschen Akademie der Wissenschaften zu Berlin, Klasse für Mathematik, 
Physik und Technik, Jahrgang 1960, Nr. 8) 


1961. 60 Seiten — 4° — DM 6,50 


Der Autor veröffentlicht Ausarbeitungen seiner Vorlesungen, die er im Herbst 1959 gehalten 
hat. Sie sollen dazu beitragen, die Schwierigkeiten der divergenten Selbstmassen und Selbst- 
‚ lJadungen in quantisierter Feldtheorie zu umgehen oder zu beseitigen. Er geht auf die Konse- 
quenzen der Lorentz-Invarianz in einer beliebigen Feldtheorie ein und charakterisiert ältere 
Versuche zur Vermeidung von Ultra-Violett-Divergenzen. Ferner bringt er ausführlich 
den wichtigen allgemeinen Beweis von Lehmann und Källen über das Auftreten von singu- 
lären Ausbreitungsfunktionen in beliebigen Feldtheorien. In zwei Kapiteln behandelt der 
Verfasser die indefinite Metrik im Hilbertraum und ihre Problematik und referiert über den 
gegenwärtigen Stand der Heisenbergschen Urmaterie-Theorie. Er geht auf die „‚axiomatische 
Methode“ in der Quantenfeldtheorie ein, die z. B. Wightman, Lehmann, Källen, Haag und 
andere stark vertreten. Abschließend versucht er die gegenwärtige Entwicklung der Quanten- 


feldtheorie kritisch einzuschätzen. 


ZEITSCHRIFT FÜR 
ANGEWANDTE MATHEMATIK UND MECHANIK 


GAMM-Tagungsheft 1960 


1960. 156 Seiten — 83 Abbildungen — 4° — DM 15,— 


In diesem GAMM-Sonderheft sind die Haupt- und Kurzvorträge enthalten, die auf der 
vorjährigen Jahrestagung der Gesellschaft für- Angewandte Mathematik und Mechanik 


in Freiberg/Sachsen gehalten wurden. 
Sie beschäftigen sich mit den Gebieten Angewandte Mathematik, Rechentechnik, Wahr- 


scheinlichkeitsrechnung und Statistik, Mechanik, Elastizitätslehre und Strömungslehre. 


Bestellungen durch eine Buchhandlung erbeten 


RER REDEKeM DE. NER LSAG .- ÜBER L.IN 


Aus 18 wissens ch: 
Wilhelm Ostwalds 
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2 Briefwechsel n mit Ludwig Boltzmann, Max P E37; L. 
| Helm und Josiah Willard Cibbe, 


'Svante Arrhenius die klassische RE EN Chemie ‚gründete 
hinaus als hervorragender Farbenforscher und als N! 

kultur- und naturphilosophischer Anschauungen srtindeh 

stand in einem ‚ausgedehnten Briefwechsel mit bedeutenden 6 

genossen seiner Epochen. Aus, ‚der umfangreichen | 

> | schers, die sich i im Besitz des Wilhelm- ÖstwaldAschirs der Deutschen AI 

Wissensehäere zu Berlin in Großbothen (Sachsen) befindet, werden in d 

ausgabe erstmalig wichtige Korrespondenzen Ostwalds- veröffentlicht, 3: 

. vorliegenden 1. Teil der Briefwechsel mit den theoretischen Physikern | ig 

mann, Max Planck und J osiah Willand Gibbs und dem Mathematiker Georg Hi 

Darin kommen vorwiegend Fragen der Thermodynamik, wie Gültigkeit und A 

barkeit der thermodynamischen Hauptsätze und darauf basierende Vorungen 

der Energetik Ostwalds zur Sprache. Die in den letzten ‚Jahrzehnten des. 

Jahrhunderts geschriebenen Briefe geben nicht nur einen interessanten Einblick in di a 

Zeit des Überganges von der klassischen zur modernen Physik, sondern vermitteln 

zugleich ein lebendiges Bild der betreffenden Bra und en 

lichen Beziehungen, die zwischen ihnen bestanden. 2 

Der mit einer Einleitung, Anmerkungen und Register versehene Band wendet sich 

hauptsächlich an Physiker, Chemiker und alle an der Geschichte der Naturwissen- 


schaften interessierten Leser. 


= 


Im zweiten, in Vorbereitung befindlichen Teil der Ausgabe werden die Korrespondenzen 
Ostwalds mit S. Arrhenius und J. H. van’t Hoff gen 


Bestellungen durch eine Buchhandlung erbeten 
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